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Resume en Anglais 


This work is devoted to the study of the existence, uniqueness, periodicity 
and positivity of solutions for two classes of delay second order differential 
equations. The first class is of neutral differential equations with time varying 
delay while the second one is of iterative differential equations with time and 
state dependent delays. The method used here is one of the most efficient 
techniques for studying these types of equations since it combines some useful 
properties of Green 's functions together with fixed point theorems to esta¬ 
blish sufficient conditions for proving the desired results. The idea of this 
technique is based on the converting of the considered equation into an inte¬ 
gral one whose solutions are recourse to a Krasnoselskii’s or Schauder 's fixed 
point theorem. In addition, by applying the principle of contraction mapping 
on the integral equation and under certain specified constraints we proved 
the existence and uniqueness of the periodic solution of the second class of 
equations. 


Keywords : Contracting principle, existence, Green's function, iterative dif¬ 
ferential equation, Krasnoselskii’s fixed point theorem, neutral differential 
equation, Schauder 's fixed point theorem, uniqueness. 
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Resume en Frangais 


Ce travail a ete devoue a l’etude de l’existence, l’unicite, la periodicite et 
la positivite des solutions de deux classes d’equations differentielles du second 
ordre a retard variable. La premiere classe est d’equations different ielles de 
type neutre avec un retard variant dans le temps tandis que la deuxieme est 
d’equations different ielles iteratives avec des retards dependants du temps et 
de l’etat. La methode utilisee ici est l’une des techniques les plus efhcaces pour 
etudier ces types d’equations, car elle combine certaines proprietes utiles des 
fonctions de Green avec des theoremes de point fixe afin d’etablir des condi¬ 
tions sufhsantes pour prouver les resultats souhaites. L’idee de cette technique 
est basee sur l’inversion de l’equation consideree en une equation integrate 
dont la solution fait recoure au theoreme de point hxe de Krasnoselskii ou 
celui de Schauder. Aussi, moyennant le principe de l’application contractante 
de Banach et sous certaines contraintes determinees, on a montre l’existence 
et l’unicite de la solution de la deuxieme classe d’equations. 


Mots cles : Principe de contraction, existence, fonction de Green, equa¬ 
tion differentielle iterative, theoreme de point fixe de Krasnoselskii , equation 
differentielle de type neutre, theoreme de point hxe de Schauder, unicite. 



L 


Resume en Arabe 
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Typographie et Acronymes 


Typographie 

r designe le retard dans l’ensemble de cette these et est, soit une fonction, 
soit une constante suivant les modeles et les exemples. 

H designe 1’argument devie dans l’equation a arguments devies. 

T designe la periode. 

Les titres des references, les noms des savants, les termes techniques ou les 
mots que nous voulons souligner sont indiques en italique. 

Les notations entre guillemets correspondent a l’emploi de termes techniques 
que nous voulons faire ressortir. 

Acronymes 

Abreviation Signification 

CSH cellule souche hematopoietique. 

dde23 solveur de Matlab, permet de resoudre les EDR a retard constant, 

ddensd solveur de Matlab, permet de resoudre des equations differentielles, 

de type neutre oil le retard est independant de I’etat. 
ddesd solveur de Matlab, permet la resolution des equations different ielles, 

a retard variable. 


IX 



Acronymes et Typography 


ddefile 

fonction au second membre de l’Squation differentielle a retard. 

ddefun 

fonction au second membre de l’Squation different ielle a retard. 

DELSOL 

solveur de Fortran, permet de rSsoudre les equations differentielles, 

a retard. 

EDF 

equation differentielle fonctionnelle. 

EDN 

equation differentielle de type neutre. 

EDO 

equation differentielle ordinaire. 

EDP 

equations aux dSrivSes partielles. 

EDPR 

equations aux derivees partielles a retard. 

EDR 

equation differentielle a retard. 

DELSOL 

solveur de Fortran, permet de resoudre les equations different ielles, 

a retard. 

Fortran 

langage de programmation compile utilise principalement pour, 

le calcul scientifique. 

history 

fonction historique. 

lags 

vecteur colonne represente les retards. 

Maple 

langage interprets de programmation utilise principalement pour, 

le calcul formel. 

Matlab 

langage de programmation interprets utilisS principalement pour, 

le calcul matriciel et numerique int-eractif 

p.p 

Presque partout. 

tspan 

intervalle d’intSgration. 

RADAR5 

solveur de Fortran, permet de rSsoudre les Squations different ielles, 

a retard. 

RETARD 

solveur de Fortran, permet de rSsoudre les Squations different ielles, 

a retard. 


x 
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Notations 


R 

M+ 

M n 

N 

N* 

[a,b] 

(a, 6] 

(a, 6) 

C = C([a,b] ,M) 

H 1 = H 1 (([—r, 0], M n )) 
L 1 ((0, +oo), M + ) 


Pt 

C 

1 

© 


Ensembles et nombres 

ensemble des nombres reels. 

ensemble des nombres reels positifs ou mils. 

espace euclidien de dimension n. 

ensemble des entiers naturels. 

ensemble des entiers naturels non nuls. 

intervalle ferine de M d’extremites a et b. 

intervalle semi-ouvert de M d’extremites a et b. 

intervalle ouvert de M d’extremites a et b. 

espace des fonctions continues de [a, b] dans M. 

espace de Sobolev d’ordre 1. 

espace de Lebesgue des fonctions sommables, 

sur (0, +oo). 

ensemble des fonctions continues et T—periodiques. 
ensemble des nombres complexes, 
ensemble des fonctions continues et T—periodiques. 
ensemble des fonctions continues et T—periodiques. 
telles qu’il existe K>0etL>0:K^ip^L. 



Notations 


Pt(L,M) 


x 


t £ M 
B' (a, r ) 

C 0 (R) 

K 


ensemble des fonctions x de Pt de norme ||x|| < L, 

telles que \x (t 2 ) — x (fi)| < M \t 2 — t x \, Vti,t 2 £ K. 

x = — est la derivee de la variable x par rapport au temps et, 
at 

c’est la derivee a droite dans les resultats d’existence et d’unicite 

dans la deuxieme section du deuxieme chapitre. 

variable temporelle. 

boule fermee de centre a et de rayon r. 

espace de toute les suites reelles avec lim Q = 0. 

i ->oo 

IK = M on C. 


Fonctions 


|-| valeur absolue. 

Il/Hoo norme de la convergence uniforme et se definit par sup \f (x)\ . 

max fonction maximum. 

min fonction minimum. 

sup fonction supremum. 

expx fonction exponentielle de 

G ( t , x) fonction de Green. 


xii 
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Introduction generale 


D ans la nature, plusieurs phenomenes sont gouvernes par une classe 
d’equations differentielles a retard. Comme leur non l’indique, ces 
equations notees en abrege EDR modelisent generalement revolution des va¬ 
riables dependant non seulement des valeurs actuelles mais dependent aussi 
irreductiblement des valeurs prises dans le passe autrement dit elles tiennent 
compte de l’effet du passe dans la prediction du futur c’est pour ga qu’elles 
sont souvent appelees "les equations a memoire ", "les equations hereditaires ", 
" equations a post.-effet", ou encore " equations a argument devie". 

Le retard est generalement une constante positive, une variable depen- 
dante continument du temps ou de l’etat ou distribue et se traduit comme 
un temps necessaire pour que le systeme reponde a une certaine evolution, 
ou parce qu’un certain seuil doit etre atteint avant que le systeme ne soit 
active. 

Sa signification differe d’un modele a un autre, il peut representer la duree 
du cycle cellulaire dans la division cellulaire, le temps de gestation ou de 
developpement dans la dynamique des populations, la periode d’incubation 
d’une maladie contagieuse en epidemiologie etc. 

L’introduction d’un tel retard a fait cauchemarder plusieurs chercheurs 
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pendant un siecle car sa presence peut induire des oscillations indesirables 
ou peut causer l’instabilite du systeme et rend la theorie difficile et conduit 
a faire appel a un arsenal mathematique tres vaste. Malgre qa. il a fait rever 
d’autres aussi comme Picard qui a souligne lors de la conference Interna¬ 
tionale des mathematiciens a Rome en 1908, l’importance de modeliser les 
systemes physiques par des equations retardees et il a affirme : 

Les equations differentielles de la mecanique classique sont telles qu’il en 
resulte que le mouvement est determine par la simple connaissance des posi¬ 
tions et des vitesses, c’est-a-dire par I’etat a un instant donne et a I’instant 
infiniment voisin. 

Les etats anterieurs n’y intervenant pas, I’heredite y est un vain mot. 

L’application de ces equations ou le passe ne se distingue pas de Vavenir, ou 
les mouvements sont de nature reversible, sont done inapplicables aux etres 
vivants. 

Nous pouvons rever d’equations fonctionnelles plus compliquees que les 
equations classiques parce qu’elles renfermeront en outre des integrates prises 
entre un temps passe tres eloigne et le temps actuel, qui apporteront la part 
de I’heredit eQ 

D’autre etudes ont montre egalement que l’introduction volontaire du 
retard peut ameliorer la stabilisation des systemes : amortissement et stabili¬ 
sation P], resonateurs retardes [45j, rejet de perturbation [46], [82], controle 
de cycle limite non lineaire [1]. 

La litterature concernant les equations fonctionnelles a retard est abondante, 
on ne peut pas citer le grand nombre d’ouvrages dedies a cette thematique 
mais le lecteur interesse pourra se referer aux livres de Hale-Cruz [40] . Hale 
0B. Hale- Verduyn Lunel m, Kolmanovskii-Myshkis [08j |49j et Smith 
1 J.K. Hale, History of delay equations, dans |J5], 1 — 28. 
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Plusieurs approches ont ete utilisees pour traiter ce type d’equations 
comme la methode des etapes, la methode de Runge Kutta, les transfor¬ 
mations de Laplace et de Schroder et la theorie de point fixe qui est un 
outil precieux pour la resolution des equations differentielles non lineaires en 
transformant le probleme d’existence en un probleme de point fixe. 

Le fil rouge de cette these est d’appliquer cette derniere technique pour 
etudier certaines classes d’equations differentielles non lineaires a retard lar- 
gement rencontrees dans la modelisation de plusieurs phenomenes dans des 
domaines tres varies comme la physique, la chimie, la biologie, l’ecologie, 
etc. Plus precisement, nous nous interesserons a appliquer une technique qui 
combine un des theoremes de point fixe de Banach, de Schauder ou de Kras- 
noselskii avec quelques proprietes utiles des fonctions de Green de certaines 
equations differentielles lineaires pour arriver aux resultats desires. 

Done, notre travail de these s’interesse aux problematiques suivantes : 
l’existence, l’unicite, la periodicite et la positivite des solutions de deux 
classes d’equations differentielles a retards variables, la premiere est de type 
neutre avec des retards variant dans le temps et la deuxieme est de type 
iterative avec des retards dependants de l’etat. 

On a structure ce manuscrit en quatre chapit.res et une annexe sur la syn- 
taxe des principaux solveurs de Matlab et Maple qui permetteut la resolution 
numerique des equations differentielles a retard. 

Dans le premier chapitre, nous avons juge utile de rappeler brievement 
quelques notions generales et certains resultats fondamentaux indispensables 
pour le reste de la these comme la fonction de Green unidimensionnelle et 
le theoreme d’ Ascoli-Arzela ou la partie principale du chapitre est dediee 
a la presentation des theoremes de point fixe, le theoreme de point fixe de 
Banach, de Schauder et celui de Krasnodelskii. 
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Le deuxieme chapitre est un etat de Part sur les equations different idles a 
retard. Apres un bref historique et une classification de ce type d’equations, 
nous passons a un rappel de quelques notions de base et certains resultats de 
la theorie d’existence et d’unicite. Ensuite, nous nous interesserons plus par- 
ticulierement a presenter quelques methodes de resolution comme la methode 
des etapes, la methode de Runge-Kutta et la transformation de Laplace. Nous 
proposerons pour conclure des moddes a retard modelisant cinq phenomenes 
dans l’ecologie, l’electrodynamique classique, la biologie et l’electronique. 

Enfin, l’essentiel de cette these est presente dans le troisieme et le qua- 
trieme chapitres qui correspondent a des articles publies [13 US] ; seules les 
references ont ete regroupees dans une bibliographie commune. La methode 
utilisee dans ces deux chapitres est tres efhcace car elle combine certaines 
proprietes utiles des fonctions de Green et la technique du point fixe. Dans 
le troisidne chapitre, on va utiliser le theoreme de Krasnoselskii pour etablir 
l’existence des solutions pd’iodiques positives de l’equation differentielle non 
lineaire du second ordre de type neutre suivante : 

^px(t)+p(t)^-x(t) + q(t)x(t) = ^pQ(t,x(t-T(t))) 

+f(t, hi (. x{t )), gi {x(t - r(t)))), 

ou p et q sont- des fonctions positives continues. La fonction Q : IxK —> M 
est differentiable, hi, g\ : M —> let/:lxlxl —* M sont des fonctions 
continues par rapport a leurs arguments. 

Le travail expose dans ce chapitre est realise en collaboration avec Pr 
Ahcene. Djoudi. et Dr Abdelouaheb. Ardjouni., et est publie dans le journal 
Nonlinear Dynamics and Systems Theory [IB]. 

Le quatrieme et dernier chapitre est dedie a l’etude de l’equation diffe- 
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rentielle iterative du second ordre suivante : 

^px(t)+p(t)^x(t) + q(t)x(t) = ^g(t,x(t),x [2 \t),...,x [n] (t)) 

+h(t,x(t),x [2] (t),...,x [n] (t)) , 

ou x (t) = t, xM(i) = x(t), x^(t) = x (x(t )), ...,x^(t) = x^ 1 } (x(t)). Les 
resultats exposes dans ce chapitre est une publication dans le journal Acta 
Mathematica Universitatis Comenianae [16]. 

On termine cette these par une conclusion generale sans doute non ex¬ 
haustive ainsi que les perspectives de nos recherches futures. 
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A cause de leur grande importance dans notre etude, ce chapitre intro- 
ductif vise a presenter quelques notions de base necessaires pour le 
developpement ulterieur de cette these et on insistera en particulier sur les 
fonctions de Green unidimensionnelles, au theoreme d 'Ascoli-Arzela et aux 
theoremes de point fixe de Banach , de Schauder et celui de Krasnoselskii. 
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Chapitre 1. Notions Preliminaires 

1.1 Fonctions de Green uni-dimensionnelles 

La fonction de Greer^j qui est une solution elementaire d’une equation 
lineaire que ce soit differentielle ou aux derivees partielles, est un outil tres 
puissant pour la resolution de ces equations et son interet principal reside 
dans la representation integrate de la solution. 

On va aborder cette notion que tres legerement ici, juste pour rappeler 
les conditions de son existence. 

Soit l’equation different ielle du n ieme ordre suivante : 

a 0 it) y {n) ( t ) + ai (t) y n 1] ( t) + ... + a n (t) y(t) — h (t), Vt G [c, d] , (1.1) 

ou les fonctions ao, a i,..., a n , h sont continues sur [c, d ], ao (t) f 0 sur [c, d) 
avec les conditions aux limites suivantes : 


oi\y (c) + a{\' (c) + ... + a± l ^y( n b (c) + fi x y (d) + ... + ^ y( n b (d) = 0 

a 2 y (c) + a^y' (c) + ... + ( c ) + p 2 y (d) + ... + /jO-b^n-b (^) = o 


a n y (c) + offy' (c) + ... + oif 1 l) y^ n 1} ( c ) + f3 n y (d) + ... + /3 ( n n 1} y (n 1} (d) = 0. 

( 1 . 2 ) 


L’equation homogene de (1.1) est : 


a 0 (t) V (n) if + oi [t) y n 11 it) + ... + a n it) y (t) = 0, Vf G [c, d] . (1.3) 

1 Cette fonction doit son nom au physicien britannique George Green (1793 — 1841). 
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Chapitre 1. Notions Preliminaires 


On suppose que le rang du tableau A ci-dessous est n : 



a i 

«! 1} ■ 

. (5, . 

.. 

A = 

a 2 

«2 1} • 

• P 2 • 

I3t l) 


(y, n 


• Pn ■ 

p { r' ] 


Theoreme 1.1 


a) Si l’equation homogene de (1.1) avec les conditions frontieres (1.2) n’ad- 


met que la solution triviale, alors il existe pour ce probleme une fonction de 
Green unique G (t, s ) on suppose o 0 (t) ^ 0, Vi G [c, d]. 

b ) on suppose realisees les conditions ci-dessus, alors Vequation non homo¬ 


gene (1.1) admet une solution unique verifiant les conditions frontieres (1.2) 
donnee par : 


v(t) = G(t,s) g(s)ds. 


(1-4) 


Supposons que le probleme aux limites (1.3)-(1.2) n’admet que la solution 
triviale y ( t ) = 0. 


Definition 1.1 


On appelle "fonction de Green " ou "fonction d’influence" 


de 1’equation (1.3) avec les conditions frontieres (1.2), la fonction G(t,s 


construite pour tout point s G ]c, d[, et jouissant des quatre proprietes sui- 
vantes : 

1/ G (t, s ) est continue et possede des derivees continues par rapport a t jus- 
qu’a l’ordre (n — 2) inclus pour t G [c, d] . 

2/ Sa (n — l) r " lf derivee par rapport a t presente au point t — s une discon¬ 
tinuity de premiere espece, le saut ayant la valeur ——, i.e. 

a 0 (s) 
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ou 


dG 

~dt 


(s + ,s) et 



s ) 


designent respectivement les limites a droite et a gauche lorsque t tend vers 

dG 

s de la fonction qui a t fait correspondre (t, s ). 

3/ Pour tout s fixe dans ]c,d\, cette fonction verifie la condition frontiere en 
c et la condition frontiere en d. 

4/ Pour tout s fixe dansjc, d\, la fonction qui a t fait correspondre G (t, s ) est 
solution de l’equation homogene dans chacun des intervalles [c, s ) et (s, d\. 


1.2 Theoreme d’Ascoli-Arzela 

Ci-dessous on rappelle le theoreme d 'Ascoli-Arzela qui est un outil classique 
et puissant pour montrer qu’une partie de l’espace des fonctions continues 
sur un compact est relativement compacte. 

Definition 1.2 Soit M un sous ensemble de E = C ([a, 6],M), M est dit 
equicontinu si et seulement si 

We >0, Wx G [a, b\, > 0, Wy e [a, b ], 

[||x - y\\ < 5] =>• [V/ G M, I/O) - f(y)\ < e]. 

Exemple 1.1 

Si k est un nombre reel positif, alors l’ensemble des fonctions lipschitziennes 
de E dans M de rapport k est equicontinu. En effet, pour e fixe, il suffit de 
prendre 

£ 

a — TV) 

M 

dans la definition 1.2. 


9 






Chapitre 1. Notions Preliminaires 


Theoreme 1.2 IUUJ Soit C([a, 6 ],M), — oo < a < b < +oo, I’espace des 
fonctions continues definies sur le compact [a, b] et a valeurs reelles muni de 
la norme 

IImII = max \u(t)\. 

a<t<b 

Une partie T de C ([a, b ], M) est relativement compacte dans C ([a, b ], M) si et 
seulement si elle est uniformement borne^\ et equicontinue. 

Remarque 1.1 Le theoreme d’ Ascoli-Arzela ne permet de caracteriser que 
les ensembles relativement compacts de C (Ei, IK), (avec E x compact et me- 
trique, K = M ou C), et non les ensembles relativement compacts de n’im- 
porte quel espace de Banach. 

Exemple 1.2 Soient k\ et /r 2 deux reels strictement positifs. Le sous-ensemble 
T des fonctions reelles continues sur [a, b], derivables sur }a,b[ qui verifient 

1 / (t)| < h et sup |/' (t)\ < k 2 , 

pour tout t G [a, b], est relativement compact dans C([a, 6 ],M). 

En effet, pour tout / G T. le theoreme des accroissements finis, prouve que 
pour tout to, t G [a, b] il existe c G ]to, t[ tel que 

1 /(*) - / (*o)| = |/' (c)| |t - tol - 

£ 

Done |/(£) — / (to) | < k 2 \t —1 0 \. Fixons t 0 G [a, b]. Soit e > 0 et 77 = —, 

k 2 

alors 

V£ G [a,b\, \t-t 0 \<r)^\f (■ t) - f (f 0 )| < e- 
est dite "uniformement bornee" si et seulement si 

\/x e [a, 6 ], {/(*),/£ J 7 } est borne. 
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Ce qui est exactement l’equicontinuite de T en t 0 . Comme nous pouvons 
prendre pour to n’importe quel point de [a, b], on en deduit que T est equi- 
continu. 

On a |/ (t)\ < k\ pour tout / e T ce qui implique que ||/’|| 00 < k\ et partant 

v/ G tF, f e B' (0, ki ), 


i.e., 


T C B' (0, k x ), 


d’ou la bornitude de T. 

Finalement, Comme T est borne et equicontinu, alors le theoreme d'Ascoli- 
Arzela assure que T est relativement compact. 


1.3 Theoremes de point fixe 

Comme l’attestent les tres nombreux travaux paraissant aujourd’hui au ni¬ 
veau international, les theoremes de point fixe sont des outils precieux et 
tres interessants en mathematiques, surtout pour la resolution des equations 
differentielles non lineaires. Pour resoudre un probleme par la technique du 
point fixe, nous avons besoin d’une application appropriee, d’un ensemble 
convenable apte pour contenir les solutions du probleme et d’un theoreme de 
point fixe qui donne certaines conditions sous lesquelles cette application ad- 
met au moins un point fixe. On va voir maintenant trois theoremes de point 
fixe, le theoreme de point fixe de Banach qui donne un critere general dans 
les espaces metriques complets, celui de Schauder qui est plus topologique et 
affirme qu’une application continue sur un convexe compact admet au moins 
un point fixe et finalement, le theoreme hybride de Krasnoselskii. 
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Definition 1.3 Soit (E, ||.|| s ) un espace de Banach , T : E —» E une ap¬ 
plication. On appelle point fixe de T tout point x G E tel que T(x) = x. Ce 
qui est equivalent a dire que l’equation T(x) — x = 0 possede une solution. 

Definition 1.4 Soit ( E , ||.|| E ) et (F, ||.|| F ) deux espaces de Banach et 
T : E —» F une application. On dit que T est lipschitzienne de rapport 
k^O sur E si 

\\Tx-Ty\\ F < k\\x-y\\ E . (1.5) 

En particulier, si 1 > k ^ 0. T est dite contraction ou application contractante 
de rapport k. 

1.3.1 Theoreme de point fixe de Banach 

En 1922, le mathematicien polonais Stefan Banach a prouve son celebre 
theoreme (connu aussi sous le nom de theoreme de l’application contractante 
ou theoreme de Banach- Picart^ qui garantit l’existence et l’unicite d’un 
point fixe d’une application contractante d’un espace metrique complet dans 
lui-meme. En outre, il est base sur un processus iteratif assurant que ce point 
fixe peut etre obtenu comme limite d’une suite iteree et qu’il est possible 
d’estimer la precision avec laquelle cette limite est atteinte. 

Theoreme 1.3 fWj Soit (E, Il-H^) un espace de Banach et T : E —> E 
une application contractante de constants k G [0,1 [. Alors il existe un point 
unique x 6 E tel que T(x) = x. 

Remarque 1.2 Si A est une application lipschitzienne (pas necessairement 
une contraction) mais l’une de ces iterees A p est une contraction, alors A a 

3 Charles Emile Picard fut le premier a utiliser le theoreme de point fixe de Banach 
dans une methode d’approximations successives. 
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un seul point fixe. 

En effet, soit x V unique point fixe de A p on a 

A p {A{x)) = A{A p (x)) =A(x), 

ce qui convient a dire que A (x) est aussi un point fixe de A p et grace a 
l’unicite A (x) = x. 

Ce resultat est valable pour tous fes types de contractions qui assurent Funi¬ 
cite du point fixe. 

1.3.2 Theoreme de point fixe de Schauder 

Le theoreme de point fixe de Schauder elabore en 1930, assure l’exis- 
tence d’au moms un point fixe pour une application continue sur un convexe 
compact dans un espace de Banach. 

Theoreme 1.4 J75|/ Soit M un sous ensemble convexe, ferme, borne et non 
vide d’un espace de Banach E et T : M —» E une application compacte. 
Alors T possede un point fixe. 

Remarque 1.3 Si M est compact et convexe, il sufRt que T soit continue 
pour avoir un point fixe pour T. 

1.3.3 Theoreme de point fixe de Krasnoselskii 

En 1954, Kranoselskii combina le theoreme de point fixe geometrique de 
Banach et le theoreme de point fixe topologique de Schauder en un theoreme 
hybride qui affirme que dans un convexe compact, toute application qui s’ecrit 
sous forme d’une somme de deux applications dont l’une est une contraction 
et l’autre est compacte admet au moins un point fixe. 
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Principe du theoreme de point fixe de Krasnoselskii [17] 

Selon Smart ra. Krasnoselskii etudia un article de Schauder m et formula 
le principe suivant; 

L ’inversion d’un operateur differentiel perturbe peut engendrer une somme 
de deux operateurs, une contraction et un operateur compact. 

Pour mieux comprendre cette observation de Krasnoselskii, on considere 
l’equation differentielle suivante : 


x'(t) = -a(t)x(t) - g(t,x(t)), 


( 1 . 6 ) 


ou a(t + T) = a{t) et g(t + T,x) = g(t, x) pour un certain T > 0. 

On peut transformer cette equation sous une autre forme en ecrivant, 

x'(t)exp (^J a(s)ds^j = —a(t)x(t)exp (^J a(.s)ds^j —g(t, x(t))exp (^J a(s)ds 

(1.7) 

Par consequent, 


x(t)exp a(s)ds\ =—g(t, x(t))exp ( / a(s)ds 


( 1 . 8 ) 


Une integration de t — T a t donne 


^ ^ x[u)exp J a(s)ds S j du = — J g{u,x{u)) ^exp j a(s)ds^j du. 

(1.9) 


Ainsi, 


(t) = x(t — T)exp a^ds'j — J g(u,x(u)) exp J a(.s)ds^ du. 

( 1 . 10 ) 


Si on suppose que 


exp ( — / a(.s)ds ) = a < 1, 
Jt-T 
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et si (E, ||.|| E ) est l’espace de Banach des fonctions p : 


continues et 


T— periodiques, alors l’equation (1.10) peut se mettre sous la forme 


<p(t) = (Bp)(t) + (Ap)(t) (:= V<p{t)), 


( 1 . 11 ) 


ou B est une contraction de constante a < 1 et A est une application com- 
pacte. Cet exemple montre bien la naissance de l’application Vp := Bp + Ap 
qui s’identifie a une somme d’une contraction et une application compacte. 


La recherche d’une solution pour (1.10) exige done un theoreme adequat qui 


s’applique a cet operateur hybride V et qui peut conclure l’existence d’un 


point hxe qui sera a son tour, solution de l’equation initiale (1.6). Krasno- 


selskii trouva la solution en combinant les deux theoremes de Banach et celui 
de Schauder en un seul theoreme hybride mais puissant qui porte son nom. 
II etablit le resultat suivant : 


Theoreme 1.5 J73j/ Soit ( E , ||.|| B ) un espace de Banach et soit M une partie 
convexe, fermee et non vide de E. On suppose que A, B : M —> E sont deux 
applications satisfaisant : 

( i) Ax + By G M, Vx, y e M; 

(ii) A est continue et A est compacte; 

(Hi) B est une contraction de constante a < 1. 

Alors, il existe x* G M tel que 

Ax* + Bx* = x*. 


Preuve. D’une part, la 3 eme condition donne 


< 


II (x-y) - (Bx-By)\\ 
\\x - y\\ + ||Bx - By\\ 


< \\x - y\\ + a\\x - y\\ 


< (l + a)\\x-y\\, (1.12) 
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et d’autre part 


\\(I-B) (x)-{I~B) (y)|| 


= \\(x - y) - (Bx - By)\\ 

> \\I (x - y)\\ - \\Bx - By\\ 

> \\x - y\\ - a\\x - y\\ 

> ( 1 — a) \\x**y\\ ■ 


(1.13) 


De (1.12) et (1.13), il resulte que 


(1 - a) ||x - y || < ||(/ - B) (x) - (/ - B) (y)|| < (1 + a) ||x - y||. (1.14) 


Cette double inegalite montre que (/ — B) : M —> (/ — B) M est continue 
et bijective. Done, (/ — B) 1 existe et elle est continue. 

Posons f/:=(/ — B)~ l A, il est clair que U est une application compacte, 
puisque U est une composition d’une application continue et une application 
compacte. 

En vertu du theoreme de Schauder , U admet un point fixe, i.e. 

3 x <G M tel que (/ — By 1 Ax = x. (1.15) 


Ceci equivaut a dire 


et la preuve est achevee. ■ 


Ax + Bx = x, 


(1.16) 


Remarque 1.4 Si A = 0, le theoreme se resume au theoreme de Banach et 
si B = 0 alors le theoreme n’est autre que le theoreme de Schauder. 

Remarque 1.5 j61| Dans une des versions generalisees de ce theoreme, la 
condition (Hi) est remplacee par la condition suivante : Il existe une fonction 
continue a droite 

0 : [0, +oo) —> [0, +oo), 
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telle que 4>(r) < <fi (s) pour 0 < r < s, (f> (V) < r si r > 0, et 

|| Bx — By || < (f) (||x — y||) pour tout x, y G M. 

Remarque 1.6 [61J Dans le cas ou a = 1 dans la condition (in) , il n’y a 
pas de point fixe. En efiet, void un contre-exemple. 

Soit l’espace de Banach C 0 (M) de toutes les suites reelles avec lim = 0; 

i -xx) 

i.e., 

Co (M) = |c = (CiW , Ci e M : .lin^Ci = o} , 

muni de la norme 

IICII = max{|Ci| : i = 1 , 2 ,...}, C = (Ci) ieN * e C 0 (M). 

On considere l’application A definie par 

A: B'( 0,1) —> C 0 (K) 

C 1 * AC, // , 

ou B'(0, 1) designe la boule unite fermee de Co (M) et 

Vi = 2 ' 1 (1 + IICII), Vi = ( 1 - 2_l_1 )Ci-i Pour i = 2, 3,... 

Puisque \rh\ < 1 et |? 7 j| < |(d| < 1, alors, on obtient 

.A: £'(0,1) — >B\ 0,1). 

En outre, si ( = (C i) ieN * et ( = sont deux elements distincts de 

B'( 0,1), on aura 

AC - A(' =rj-ri = - r/A 

V /ieN* 

OU 

Vi ~ Vi = 2' 1 (i + IICII) - + lie'll) = s^dlCII - IIC'II), 
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et 

Vi-Vi = (l-2“^ 1 )Ci_i-(l-2^“ 1 )Cl_i = (1—2 - * -1 ) (Ci-i - d-i) P°ur i = 2, 3,. 
Alors 




= max 


Vi ~ Vi 


: i — 1, 2,.. 


= max 12 ^HCII - ||C ||), max <j (1 

< iic-c'ii- 




) Ci— 1 Ci—1 


: i — 2,3, 


S’il existe ( £ B'(0, 1) tel que AC, = ( alors, Ci = 2 X (1 + ||C||) > 0. De plus 
pour tout i > 2, 


ICd = |(l-2- i - 1 )C i _ 1 | = |(l-2- i - 1 )(l-2-%_ 5 


i—2 


IK 1 - 2 


— i— l-\-k 


)Ci 


k =0 

i +1 

i - 

3 =3 


> 


i—2 


‘-E 2 


—i—l+k 


k =0 


ICil 


ICil > 4IC1I, 


ce qui contredit le fait que tend vers zero quand i tend vers l’infini et 
prouve que A n’a pas de point fixe. 


1.4 Recapitulatif 

Dans ce chapitre, nous avons aborde brievement quelques rappels et no¬ 
tions de base en vue de les utiliser dans les chapitres suivants comme la 
fonction de Green et le theoreme d’Ascoli-Arzela qui donne les conditions 
suffisantes pour qu’une partie d’un espace de Banach soit relativement com- 
pacte sans oublier de presenter les trois theoremes captivants de point fixe, le 
theoreme de Banach, de Schauder et celui de Krasnoselskii qu’ils ont initie de 
nombreuses etudes et ont ete etendu dans plusieurs directions en modifiant 
certaines de leurs hypotheses. 
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L es objectifs de ce chapitre sont multiples. Apres un apergu historique 
sur les EDR, la deuxieme section sera consacre entierement a la classi¬ 
fication de ces equations retardees selon les retards cites dans la litterature 
avant de nous interesser a resoudre quelques EDR par trois methodes : la me- 
thode des etapes et de Runge-Kutta et la transformation de Laplace. Enfin, 
nous avons trouve judicieux de presenter certains modeles retardees. 
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2.1 Contexte historique 

Les equations fonctionnelles a retard constituent un champ d’etude tres 
important pour modeliser des phenomenes d’heredite rencontres en physique, 
biologie, chimie, economie, ecologie, etc. Malgre que dans la plupart des mo- 
deles, le retard est estime non signicatif et ignore pour simplifier l’etude, il a 
ete prouve que dans de nombreux cas. le retard joue un role dominant dans 
plusieurs domaines et que les modeles avec retard fournissent des resultats 
plus precis et realistes que leurs homologues sans retard. 

A notre connaissance l’apparition de ces equations remonte au I8 eme 
siecle, elle est due a J. Bernoulli, L. Euler , J.L. Lagrange , P. Laplace , S. 
Poisson et d’autres. Par exemple, en 1728 dans ses experiences sur la corde 
vibrante et en partant d’une equation aux derivees partielles de type hyper- 
bolique, Bernoulli |1IJ a trouve l’equation a retard suivante : 

y' = V (t - 1), 

mais malheureusement il decida la tenir pour fausse et il est dit qu’il y a eu 
plusieurs erreurs a deduire l’equation. 

Ce type d’equations a reste terra incognita jusqu’au debut du 20 siecle 
et les travaux pionniers qui etablissent le debut de la theorie ont ete dans 
la geometrie et la theorie des nombres et les premiers papiers traitant les 
equations fonctionnelles retardees lineaires sont dus a Polossuchin (1910) 
[FIT] et Schmidt (1911) [75] . 

Dans Page d’or de l’ecologie theorique, l’etude de ce type d’equations 
connut un essor considerable avec les series de travaux de V. Volterra [sa 
i85] sur les modeles predateur-proie et les modeles de viscoelasticite. Il a 
utilise la methode d’energie pour etudier une classe generale d’equations a 
retard non lineaires et il a ecrit dans son oeuvre majeure [M] sur le role des 
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effets hereditaires sur les modeles de la dynamique de plusieurs especes en 
interaction. 

II y a eu d’intensives recherches sur le sujet depuis 1940 (surtout en l’ex- 
Union sovietique). La regulation sur base de modeles lineaires et stationnaires 
avec retard fut abordee en 1941 par Y. Zypkin. En 1949, A. D. Myshkis 
a pose les bases de la theorie moderne des EDR. En particulier, il fut le 
premier a formuler Penonce du probleme de Cauchy pour des equations a 
retard arbitraire. 

Les annees cinquante, ont vu une explosion de la theorie qui a ete large- 
ment developpee et les EDR fait partie du vocabulaire des chercheurs tra- 
vaillant sur la viscoelasticite, les problemes mecaniques, les reacteurs nu- 
cleaires, le flux de chaleur, les reseaux de neurones, la combustion, Pinterac¬ 
tion des especes, les modeles microbiologiques, epidemiologiques ou physio- 
logiques, ainsi que beaucoup d’autres (voir [49]). 

La litterature concernant les EDR dans cette periode est abondante et de 
nombreux travaux ont etabli des resultats generiques dont on cite les resultats 
de Krasovskii iodl 53j , qui a etendu la deuxieme methode de Liapunov aux 
equations retardees, sans oublier aussi les travaux de Myshkis p2J, Krasovskii 
[531153!, Bellman et Cooke [T2] , Halanay [39] . 

Les annees suivantes ont donne naissance a un grand nombre de travaux 
dans cette direction, et surtout ceux qui concerne Panalyse de la stabilite 
des equations differentielles, avec un argument retarde comme El’sgol’ts et 
Norkin [33], Hale m, Hale et Lunel 021. Diekmann, Van Gils Lunel et 
Walther [22] , 
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2.2 Classification 

Le but de cette section est d’esquisser un tour d’horizon des principaux 
types d’equations differentielles a retard. Generalement, la pratique de la mo- 
delisation prouve que seules les equations retardees ou de type neutre sont uti- 
lisees pour representer des phenomenes reelles mais il existe un nombre consi¬ 
derable de phenomenes modelisees par des equations avancees, par exemple 
Interpretation physique de l’equation differentielle avancee suivante : 

x ( t ) = x" (t + 1), 

est celle d’une force liee par exemple a l’elongation d’un ressort et qui n’agirait 
qu’apres un certain temps. 

On a entame cette section par la presentation de la classe mere dite des 
equations different ielles a argument devie avant d’exposer les differents types 
d’equations different ielles a retard. 

2.2.1 Equations different ielles a argument devie 

Les equations different ielles a argument devie appartiennent a la classe des 
equations different ielles fonctionnelles (EDF). Elies decrivent revolution des 
variables dependent des valeurs prises dans le passe ou le futur. On distingue 
trois categories principales, les equations different ielles avancees , les equa¬ 
tions different ielles retardees et les equations different ielles mixtes. 

Equations differentielles avancees 

Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 
x ' {t) = f(t, x (t) ,x(t + n)) , /i > 0, 
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ou / est donnee. Autrement dit, la valeur de la derivee a un instant t de la 
variable liee x. ne depend pas seulement de la valeur de a: a l’instant t. mais 
aussi du futur. 

Equations differentielles a retard 

Sous leur forme la plus simple, ces EDR s’ecrivent comme suit : 

x ' if) = f (t, x (t) ,x(t-r)), t > 0, 

ou / est donnee. 

Autrement dit, la valeur de la derivee a un instant t de la variable liee x. ne 
depend pas seulement de la valeur de x a l’instant t, mais aussi des valeurs 
prises avant l’instant t. 

Equations differentielles a argument devie de type mixte (avancees- 
retardees) 

Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 

x> (t ) = / (t, x (t) , x (t + ft) , x (t - t)) , H, t > 0, 

ou / est donnee. Autrement dit, la valeur de la derivee a un instant t de la 
variable liee x. depend a la fois du passe et du futur. 

2.2.2 Equations differentielles a retard 

Les equations fonctionnelles a retards peuvent etre classees comme lineaires 
ou non lineaires, autonomes ou non autonomes, periodiques ou non ou encore 
selon les types des retards. Dans ce paragraphe, on s’interesse a donner une 
classification des EDR selon les types de retards cites dans la litterature ou 
on distingue deux classes principales, la premiere est dite "d ’equations diffe¬ 
rentielles retardees" et l’autre "d : equations differentielles de type neutre". 
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Equations differentielles retardees 

Dans cette premiere classe, on trouve les sous categories suivantes : 

Equations differentielles a retard constant Sous leur forme la plus 
simple, les equations differentielles a retard constant, discret ou ponctuel 
s’ecrivent comme suit : 


x' (t) = f(t,x(t),x(t-r)), 

ou / est donnee et r est une constante positive. A titre d’exemple, on trouve 
ce type d’equations dans le modele des mouches de Nicholson [39] . 

Equations differentielles a retard variable Le retard dans ce cas va¬ 
riant dans la variable temporelle ou bien dependant de l’etat. A titre d’exemples, 
on trouve les equations differentielles a retard depend du temps dans les mo- 
deles de transport [39] et celle qui depend de l’etat dans le modele decrivant 
revolution d’une population de poissons dont les larves consomment une 
nourriture [TO] ou le probleme de deux corps en electrodynamique classique 

[23 EH [30]. 

Equations a retard variant dans le temps 

Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 

x ' (t) = f(t,x(t),x(t-r(t))), 


ou / est donnee. 

Equations a retard variable dependant de l’etat 

Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 

x' (t) = f (t,x(t) ,x(t - t (x (t )))), 
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ou / est donnee. 

Des fois cette classification n’est pas sufRsante dans l’etude, ce qui conduit a 
ajouter des contraintes supplementaires relatives au retard ou a sa derivee et 
qui mene ensuite a definir de nouvelles sous categories d’equations a retard 
variable comme : 

Equations a retard variable arbitraire : Le retard et sa derivee ne sont 
pas limites. 

Equations a retard majore : Cette sous-categorie suppose la connaissance 
d’une valeur maximale sur le retard 

0 <r(t) < r max . 

Si t (t) = t est constant, il reste en pratique incertain et la contrainte ci- 
dessus assure un intervalle borne. Ce type de retard a ete tres largement 
aborde dans la litterature. 

Equations a retard (bi-borne) : Cette sous-categorie est moins abordee 
que le cas precedent dont on suppose que le retard verifie la contrainte 

^min C ' (t ) C T max • 

Equations a retard variant lentement dans le temps 

r (t) est une fonction derivable presque partout telle que 

t (t) < A < 1, 

qui indique alors une limitation sur la vitesse de variation du retard et que 
ce dernier varie lentement dans le temps autrement dit que les informations 
retardees arrivent dans l’ordre chronologique. 

Equations a retard variant moderement dans le temps 
r (t) est une fonction derivable presque partout telle que 

t (t) < A avec A > 1. 
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Equations a retard variant rapidement dans le temps : Dans cette 
sous-categorie, il n’y a aucune contrainte sur le retard et sa derivee. 

Equations differentielles a retard distribue Ce type d’equations a ete 
traite dans la litterature, pour la conception d’observateurs pour des systemes 
non lineaires; pour prouver la stability des systemes lineaires de type neutre; 
la stabilisation robuste des systemes neutres incertains et la commande ro- 
buste. Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 



Par exemple, la chauve souris, a la chasse, etant aveugle, elle emet des sons, 
pour utiliser les parois des grottes, afin de localise!' sa proie. L’echo obtenu 
par le rebondissement de ces cris represente le retard qui depend de l’etat, 
qui est le predateur. On trouve ce type d’equations Aussi dans le modele de 


sida [63] ou le modele de dynamique des populations presente par Volterra 


en 1934 ou il a utilise un terme de retard distribue pour examiner un effet 
cumulatif sur le taux de mortalite d’une espece. 

Equations differentielles a retard inconnu Dans ce cas, aucune hypo- 


these sur le retard n’est consideree qu’il soit constant, variable ou distribue. 


Equations differentielles a retard de type neutre 

Ces equations notees en abrege EDN differencient des EDR par le fait 


que la derivee de l’etat au temps actuel depend non seulement des valeurs 
de l’etat passe mais aussi de la derivee d’ordre le plus eleve intervenant dans 
l’equation du temps passe. 
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Sous leur forme la plus simple, ces equations s’ecrivent comme suit : 

J t [Dx ( t-T 0))] = f (t, X (t) , X (t - T (t))) , 

ou / est donnee et D est un operateur. Par exemple, on trouve ce type 
d’equations dans le modele de reseaux distribues [49:. 

2.3 Resultats d’existence et d’unicite 

Dans cette section, on presente quelques resultats classiques sur l’exis- 
tence et Funicite. 

2.3.1 Cas des equations fonctionnelles a retard 

Etant donne un nombre r > 0, C([a, 6], M n ), l’espace de Banach des fonctions 
continues definies sur [a, b] a valeurs dans M n muni de la norme de la conver¬ 
gence uniforme. Si [a, b] = [—r, 0], on pose C = C{[— r, 0], R n ) et on designe 
la norme d’un element ip E C par 


Ml = su p \<p(t)\- 

-T<t< 0 

Si t 0 G A > 0 et x E C ([£ 0 — r, to + A] , M n ), alors pour t E [t 0 , t 0 + A], on 
definit x t E C par 

x t (s) = x(t + s ), 

pour tout s E [—r, 0]. 

Definition 2.1 [4Ij Si D est un sous ensemble de I x C, / : 0 —-> M' 1 est 
une fonotion donnee et " " represente ici la derivee a droite, Fequation 

Ht)'■= f(t> x t), (2- 1 ) 
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ou 

x t {s) := x(t + s), s E [-r, 0], (2.2) 

est lin e equation differentielle fonctionnelle a retard sur D notee en abrege 
EDR et le nombre r est appele le retard. 

II est clair que le cas r = 0 correspond au cas des equations differentielles 
ordinaires. II est evident qu’une condition initiale appropriee au temps t = t 0 
exige la determination de la fonction x sur tout l’intervalle [t 0 — r, t 0 ], i.e. 

x(t) = ip(t), t G [t 0 - T,t 0 ], (2.3) 


ou [to- r, to] 


est une fonction donnee supposee continue appelee 


la condition initiale de l’equation a retard (2.1). Ainsi, l’equation (2.1) peut 
se mettre sous la forme : 


I x'(t) := f{t,x t ), t>t 0 
\ x(t) =^(t), te[t 0 -T,t 0 \, 

ou ifj est une fonction donnee continue sur l’intervalle [to — r, to]. 


(2.4) 


Definition 2.2 jH] Etant donne ip e C et to £ K, une solution de l’equation 


(2.1) est une fonction notee x(t) telle que x(t) = ^(t) si t e [t Q — r;t 0 ] et 


satisfaisant (2.1) si t G [t 0 ,t 0 + A] avec A > 0. Une telle fonction x(t) est dite 


solution de (2.1) a travers (to, w) et elle est notee souvent par 


x{t) = X (to, ip, f). 


(2.5) 
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Lemme 2.1 JJTJ Etant donnees une fonction ip G C, t 0 G M et f (t, ip) une 
fonction continue. La recherche d’une solution de Vequation (2.1\ ) a travers 
(to, ip) est. equivalente a la resolution de Vequation integrate 

xt 0 = 

x(t) = -0(0) + f tQ f(u,x u )du , t > t 0 . 


Theoreme 2.1 (Existence) J Jl\/ Pour Vequation (2.1), supposons que est 
un sous ensemble ouvert deM. x C et f G C(f2,M n ) une application continue 
sur Si (t 0 ,0) ^ Vt, alors il existe une solution de Vequation (2.1) passant 
par (t 0 , ip). 


Definition 2.3 (41J On dit que la fonction f(t, ip) est lipschitzienne par rap¬ 
port a p sur un compact K de 1 x C s’il existe une constante k > 0 telle 
que, pour tout (t, ip() G K, i — 1, 2 , on a 


\f(t, 0i) - f(t, ip 2 )\ < k \ip l ~'ip 2 \ 


( 2 . 6 ) 


Theoreme 2.2 (Unicite) JJTJ Supposons que 12 est un sous ensemble ouvert 
de K x C, f : 12 —* M n est continue et f(t,ip) est lipschitzienne par rapport 
a 0 sur tout sous ensemble compact de 12. Si (to, ip) G 12, alors il existe une 


solution unique de Vequation (2.1) passant par (to, ip). 


2.3.2 Cas des equations fonctionnelles de type neutre 

Definition 2.4 (42j Supposons que 12 est un sous ensemble ouvert delxC 
d’elements (t,ip). Une fonction D : 12 —> M n [] est dite atomique au point (3 
de 12 si D est continue ainsi que sa premiere et seconde derivee au sens de 
Frechet par rapport a ip et D,p, sa derivee par rapport a ip, est atomique en 
(3 de 12. 

l D n’est pas necessairement lineaire. 
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Definition 2.5 |42j Supposons que 12 est un sous ensemble ouvert de RxC. 
D : Q —> W l et / : 12 —> sont deux fonctions donnees continues avec D 

atomique en zero. La relation 

^D(t,x t ) = f(t,x t ), (2.7) 


est dite une equation differentielle de type neutre notee en abrege EDN. 


Theoreme 2.3 


(Existence) Si 12 est un sous-ensemble ouvert deMxC 


et (to, ip) G 12, alors il existe une solution de Vequation (2.1) passant par 
(to A) ■ 


Theoreme 2.4 JJI y (Existence et unicite) Si Q est un sous-ensemble ouvert 
delxC et f (t, ip) est lipschitzienne par rapport a tp sur tout sous-ensemble 
compact de 0, alors pour tout (to, ip) G 12, il existe une solution unique pour 


l’equation (2.1) passant par (to, ip) 


2.4 Quelques methodes de resolution des EDR 

2.4.1 Methode des etapes 

La methode des etapes (dite aussi "methode pas a pas”, "methode des 
pas" ou "methode des integrations successives") permet de resoudre nume- 
riquement les EDR et les EDN et permet par la meme occasion d’etablir 
l’existence et l’unicite de la solution. Elle fut presentee en 1965, par R. Bell¬ 
man pour des retards constants. D’autres comme El’sgol’ts et Norkin (1973) 
ont montre qu’elle restait aussi valable pour les retards variables, a condition 
que le retard ne s’annule jamais. 
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(2.8) 


Pour fixer les idees, considerons Fequation fonctionnelle lineaire retardee avec 
coefficients variables suivante : 

x' (t) = ai (t) x ( t ) + a 2 (t) x (t — t) pour tout f G [0, r] 
x (t) = ip ( t ) pour tout t G [—r, 0] . 

Dans le cas ou a\ et a 2 sont deux constantes reelles, alors Fequation est dite 
de type " Frisch-Holme" [55] . Maintenant, on va resoudre cette equation par 
la methode des etapes. Le principe de cette methode est de chercher des 
solutions sur des intervalles de type [hr, (k + 1) r] ou k G N en suivant les 
etapes suivantes : 

l ere etape : Dans l’intervalle [— r, 0], la fonction x (t) est la fonction donnee 
(p(t), done Fequation est resolue dans l’intervalle [—r, 0] et on designe par 
ay (t) cette solution. II faut signaler ici que si t G [0, r], alors t — r residera 
dans [—r, 0]. 

2 eme etape : Dans l’intervalle [0,r], si t G [0, r], alors t — r residera dans 


(2.9) 


[—r, 0], done x (t — r) = xo (t — r) dans l’intervalle [0, r] et le systeme (2.8) 
devient : 

x' ('t ) = oq (t) x ('t ) + a 2 (t) Xq (t — r) pour tout t G [0, r] 
x{0) = (p (0), 

qui est un probleme a valeur initiale pour une equation differentielle ordinaire 
(EDO) ou x 0 (t — r) = ip (t — r) est connue. Ainsi, on resout cette EDO dans 
[0, r] en utilisant la condition initiale x (0) = <p (0) et on designe par X\ (t) 
cette solution dans [0, r]. 

3 eme etape : Dans Fintervalle [r, 2 r] , le systeme devient : 


x' (■ t ) = ai (t) x (■ t ) + o 2 (t) x\ (t — t ) pour tout t G [r, 2r] 


( 2 . 10 ) 


x[t) = x x (r). 

Cette EDO avec la condition initiale x\ (r) est a son tour peut etre resolue 
pour trouver la solution x 2 (t) G [r, 2r] et ainsi de suite. 
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Exemple 2.1 Considerons le cas particulier suivant : 

dx 

— = ax{t — r), pour t £ [0, 2r] 

LLL 

x(t) — (p (t) — 1 —r <t< 0. 

l ere etape : Dans [—r, 0] 

x ( t ) = 1. 

2 eme etape : Integration dans [0, r] 

L’integration des deux membres de L’EDR de 0 a t. donne 

x (t) = a x(s — r)ds + x (0). 

Jo 

Comme 0 < s < r, alors —r < s — r < 0. Sachant que x(t) = 1 pour 
t e [—t, 0[, alors 

x(s — r) = 1, 

pour s £ [0, r], ce qui conduit a 

x (■ t ) = at + 1, 


dans l’intervalle [0,r]. 

3 eme etape : Integration dans [r, 2r] 

L’integration des deux membres de r a t, donne 

x (t) — a x(s — r)ds + x (r). 


Comme r < s < 2r, alors 0 < s — r < r. Sachant que x(f) = at + 1 pour 
t £ [0,r[, alors 

x(s — r) = a(s — r) + 1, 


pour s £ [r, 2r], ce qui conduit a 


x (t) = a 2 ^— + (a — cr 2 r) t — a 2 ^— + cr 2 r 2 — or, 
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dans l’intervalle [t,2t\. 
Finalement, on a obtenu 


t 

EDO avec condition initiate 

[0,r] 

1 

x' (t) = cr 

x(0) = 1 

[r, 2r] 

1 

x' (t) = a 2 (t — t) + a 

x (r) = ar + 1, 


et la solution est donnee par 


x{t) = 


Si on prend a = 


at +1 0 < t < r 

t 2 r 2 

a 2 — + (a — a 2 r) t — a 2 — + a 2 r 2 — ar r <t < 2r. 

— 1, r = 1 avec la meme condition init.iale tp(t) — 1, on aura 


x(t ) 


1-t 0 < t< 1 

+ 2 3 

-2f + - 1 < t < 2. 

2 2 “ “ 



Fig. 2.1 - Representation graphique de l’equation x' (t) = —x(t — 1) avec 
(p(t) — 1 (plein) et tp (t) — — 1 (pointille) [69]. 


33 























Chapitre 2. Equations differentialles a retard 


Exemple 2.2 On considere le systeme suivant : 

x(t) = Ax(t ) + Bx(t — r) + Mx(t — r), p.p. t > 0, r > 0, (2.11) 

avec la condition initiale 


x 0 (9) = 0(0), 0 £ H 1 = H 1 ([-r, 0] ,M n ), - r < 9 < 0. (2.12) 


l er cas : Si A = 0, B et M sont des matrices constantes de dimension n x n. 


la solution de (2.11) est 


x 


(' t) = x(0) + f*[Bx(s — t) + Mx(s — r)]ds, t > 0, (2-13) 


mais dans l’intervalle t £ [0, r] 


x(. -r) = 0(.). 


Le changement de variable 9 = s — r et la condition (2.12) sur l’integrale 


(2.13), nous donnent 


x(t) = z(0) + S‘- r \B<t>(e) + M<t>(9)\d9. 


(2.14) 


Alors, l’utilisation de (2.13) de nouveau, on peut calculer x (t) dans l’inter- 


valle t £ [r, 2r] et ainsi de suite. 


2 eme cas : Si yl 0 0 dans (2.11), la solution est plus facile a calculer si on 


reecrit (2.11) sous la forme 


— [x(t) — Mx(t — r)] = Ax(t) + Bx(t — r), t > 0, r > 0, (2-15) 

Uib 

avec xt £ C ([—r, 0], M n ) et la condition initiale xq (6) = 0(0), 0 £ C ([—r, 0], M n ), 
—r < 6 < 0. Alors il existe une solution unique x continue dans (—r, oc) don- 
nee par 

x{t) = Mx(t— r)+e j4t [0(O) — C0(—r)] + f 0 e A( ' t ~ s \AM+B)x{s — r)ds. (2.16) 
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On note aussi que la solution generale de (2.11) —(2.12) est pour t > 0 


x(t) = e At (f)( 0) + fge A ^ s ^[Bx(s — r) + Mx(s — r)]ds. 


(2.17) 


Ceci peut etre verifie par derivation directe de (2.16) et (2.17). 


2.4.2 Methode de Runge Kutta 


Les methodes de Runge Kutta font partie des methodes d’analyse numerique 
a un pas utilisee pour la resolution des equations different idles ordinaires. 
Une methode de Runge Kutta de rang q est une methode a un pas avec 

q 

$(ti,h,Xi) = y^'Wjkj 
3 =1 


kj(ti) 


f 




l<j<q{2.18) 


ou Wj , Cj, aji sont des constantes. 

Elle est 

explicite si = 0, Vi > j 

semi-explicite si = 0, Vi > j 
implicite dans les autres cas. 

Puis on a la relation suivante, entre le rang q et l’ordre p de la methode 


q ^ 4 q = p 
q > 4 p ^ q. 

Par exemple 

q = 1 on a la methode d’ Euler, 

q = 2 on a la methode d’ Euler amdioree par exemple, 
q = 4 on a la methode " classique " de Runge Kutta, 
q = 5 on a la methode de Runge Kutta-Merson par exemple. 

La methode que Ton a choisie est la methode" classique" de Runge Kutta de 
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rang 4 et d’ordre 4 : 


x i+l ~ x i + X (M + 2/^2 + 2/^3 + £ 4 ) , 
0 


(2.19) 


avec 


h 

k 2 

h 


f(U, x i ) 


, ( h 
J Mi + 2’ Xi + 

, ( h 
J Mi +2’ Xi + 



k A = f(ti + h, Xi + hk 3 ). 


( 2 . 20 ) 


Adaptation au probleme avec retard La methode "classiqne" de Runge 
Kutta est extensible au cas des equations differentielles hereditaires [68]. En 
effet, considerons l’EDR avec condition initiale suivante : 


x(t) = f ( t, x(t),x(t — r)) t ^ 0, 
x(t) = 9(t) — t ^ t ^ 0. 


Prenons h tel que r 


mh, m G N*. On a alors (2.19) avec 


( 2 . 21 ) 


k\ f (tii x ii x i—m) 1 


k 2 

h 

k A 


f 



h 

2 


h, x 1 

i“ m+ 2 


? 


„ | h h 

J \U + -,Xi + -k 2l x 


‘- m+ 2. 


f (f j + h, x t + . 


? 


Les valeurs intermediaries necessaries sont calculees a l’aide d’une interpola¬ 
tion polynomiale de la solution entre les noeuds i — m et i — m + 1. 


Remarque 2.1 On pent resoudre numeriquement, les equations differen- 
tielles a retard en adaptant par exemple les solveurs : 
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DELSOL, RETARD et RADAR5 bases sur la methode de Runge-Kutta 
(Fortran), dde23, ddesd et ddensd bases sur la methode de Runge-Kutta 
(Matlab) et la fonction Lambert W base sur la methode des etapes (Maple). 


2.4.3 Methode de Laplace 


On va resoudre l’equation de type Frisch-Holme suivante : 

y' (t) — c y(t) + b y (t — r), pour t > 0. 
y(t)=h(t) Vte[-r,0] 


( 2 . 22 ) 


Estimation exponentielle 


On va estimer la solution du probleme (2.22) par le theoreme suivant 


Theoreme 2.5 Soit y(t) une solution de I’equation (2.22). Alors il 
existe deux constantes positives a et (3 telles que 


\y(t)\ < a||/i||exp(/3t), 


(2.23) 


ou ||h|| = sup |h (0)|. 

—t< e»<o 

Preuve. On integre l’equation de type Frisch-Holme (2.22) de 0 a t et on 
arrive a 


y'(s)ds = c y(s)ds + b y(s — r) ds, 


i o 


'0 


'0 


done 


f*t—T 


y(t) = h( 0) + cj y(s)ds + bJ y(s) ds, 

ce qui implique que 

\y{t)\<\\h\\ + \c\ I \y(s)\ds+\b\ [ |y(s)| ds, 
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\y(t)\ < 


< 



En utilisant le lemme de Gronwall, on obtient 


\y(t )I < (1 + \b\ t) ||h|| exp ((|c| + \b\)t), 


< a ||/i|| exp (fit), 


avec a — 1 + \b\ r, fi = |c| + \b\ . m 

L’equation caracteristique 

L’equation caracteristique d’une equation differentielle homogene lineaire a 
coefficients constants est obtenue a partir des solutions non triviales de la 
forme exp (At). L’equation 


y'(t) = cy(t) + by(t- t), 


possede une solution non triviale de la forme exp (At) si est seulement si 


Aexp(Af) = cexp(At) + 6exp(A(f — r)). 


Ceci implique que 


A = c + 6exp(—rA), 


ou encore 


/(A) = A — c — 6exp(—rA) = 0. 


(2.24) 


/(A) = 0 est appelee l’equation caracteristique de (2.22). 
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La solution fondamentale 


Soit Y(t) une solution de l’equation (2.22) ayant pour condition initiate 


h(t) = 


0 si t < 0, 
1 si t — 0. 


(2.25) 


Y(t) est appelee solution fondamentale du probleme (2.22). 


Theoreme 2.6 La solution fondamentale Y(t) du probleme (2.22) est 
donnee par sa transformee de Laplace 

1 


C(Y{t))(\) = 


/(A)' 


(2.26) 


Aussi pour 7 > (3 (/3 = |c| + \b\) 


Y(t) = 


1 f exp(At) 


'(7) 


/(A) 


-d\ = lim -— : 


t — >oo2Tri J c _ iT /(A) 


dX. 


(2.27) 


Preuve. Puisque Y(t) satisfait l’estimation exponentielle (2.23) alors £(Y) 
existe {C(Y) est la transformee de Laplace de Y). On a 

POO POO POO 

/ exp (—A t)y'(t)dt= / exp (—A t)cy(t)dt+ / exp (—A t)by(t — r)dt. 
Jo Jo Jo 

^ -V-' 

J 

(2.28) 

On integre J par parties pour avoir 

POO 

J = [exp (—At) y(t)](f + / A exp (—At) y(t)dt. 

Jo 

On remplace par Y(t), on obtient 

POO 

J = [exp (—At) F(f)]“ + / A exp (—At) Y(t)dt, 

Jo 

POO 

= —1+1 A exp (—At) Y[t)dt. 

Jo 
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L’equation (2.22) devient 


-1 + / Aexp(—A t)Y(t)dt = 


/ exp (—A t)cY(t)dt + 

Jo 

poo 

+ / exp (—At) cY(t — r)dt. 
Jo 


Done, 


-1 + \C (Y (t)) (A) = cC{Y(t)) + b exp (—A(t + r))Y(t)dt 


= cC (Y (t)) (A) + frexp (— Ar) x 


+ / exp (—A t)Y(t)dt 


exp (—At) Y(t)dt 


i.e., 


On conclut que 


= c£ (Y (t)) (A) + b exp (-Ar) C (Y (t)) (A) + 
+6exp(—Ar) J exp (—At) Y(t)dt, 

= cjC (Y (t)) (A) + b exp (-Ar) C (Y (t)) (A). 


(A — c — feexp (—Ar)) C (Y (t)) (A) = 1. 


£(F(t))(A) = 


A — c — frexp (—Ar) /(A) 


La preuve est achevee. 


Representation integrale de la solution 


On va representer la solution y(t) du probleme (2.22) ayant pour condition 
initiale h(t) en fonction de la solution fondamentale Y(t). 


Theoreme 2.7 12&J La solution y(t ) du probleme (2.22) peut etre represen¬ 
tee sous la forme 

y(t) = Y (t) h(0) + b[ Y(t-9-T) h{9)dd. 
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Preuve. On applique la transformation de Laplace a chaque terme du pro- 


bleme (2.22), on obtient 

poo poo poo 

I exp (—A t)y\t)dt — cj exp (—At) y(t)dt + b exp (—At) y(t — r)dt. 

_ Jo Jo 

(2.29) 


I 


On integre / par parties et on trouve 


/ = [exp(-At)y(t)]“ + / A exp (-At) y(t)dt, 

Jo 

poo 

= — h(0) + / A exp (—At) y(t)dt, 

Jo 

= ~h(0) + XC (y) (A). 


L’equation (2.29) devient 

poo 

—h(0) + XC (y) (A) = cC (y) (X) + b exp(-A (t + t)) y(t)dt, 

Jo 

= cC (y) (A) + 6exp(—Ar) 


x 


exp (—At) y(t)dt + / exp (—At) y(t)dt 


= cC (y) (A) + 6exp(-Ar)£ (y) (A) 

[0 

+6exp(—Ar) / exp (—A t)y{t)dt. 

Done, 

h (A) £ (r/) (A) = h (0) + 6exp(—Ar) J exp (—At) h(t)dt, 
ceci implique que 


£ (y)(A) = 

Done, 

y(t) = 


/(A) L 

exp (At) 


h , (0) + 6exp(—Ar) x / exp (—A 9)h(6)d6 


(7) 


/(A) 


/?. (0) + 6exp(—Ar) x / exp (—A0) h(6)d6 


dX, 


■ '««+//?g?i 


exp(—Ar) x / exp (-X6) h{9)d6 


dX. 


h 
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Pour calculer I\. On definit une fonction uj : [—r; +oo[—» [0,1] telle que 


co(d) = 


0 si 9 > 0, 

1 si 9 < 0. 


r»+oo 


Ii = j exp(—Ar) exp (-X6) h(9)u {9) d9, 
En posant s = t + 9, on obtient 

poo 

1 1=1 exp(—A s)h(s — t)u ( s — r) ds, 
= £ (h(-T +-)u (-T + ■)) ( A). 


Maintenant, on va calculer 


[ 6^/irfA. 

'( 7 ) J 


f b -TTjvr lidX = b f exp(At)£ (Y) (A) C (h(-r + -)u (—r + •)) (A) d\, 

2 (7) h J l A J 2( 7 ) 

= bf exp(Af)£ (Y * [h(—r + -)uj (—r + •)]) (A) dX, 

2 ( 7 ) 

= b (Y * [h(-r + -)cj (-t + •)]) (t), 

= b Y (t — s ) [h(—r + s)o; (—r + s)] ds, 

Jo 

= b Y (t — s) h(—r + s)ds, 

Jo 

car uj (—r + s) = 0 dans [r, t]. Finalement {9 = s — r) 

y(t) = Y (t) h{ 0) + bj Y (t- 9-r) h(9)d9. 

Ce qui acheve la demonstration. ■ 
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2.5 Modeles d’equations differentielles a re¬ 
tard 

2.5.1 Equation de Nicholson (mouches-moutons) 

En 1950, le celebre biologiste entomologiste australien Alexander J. Ni¬ 
cholson a mene une longue serie d’experiences visant a en apprendre davan- 
tage sur des populations de mouches devoreuses de viandes responsables de 
90% des myiase^] ovines qui menacent les elevages de plusieurs pays comme 
VAustralie , la Nouvelle Zelande et VAfrique du sud [651 [61. [56]. 

Cette mouche diptere ayant un corps rond a ovale de longueur varie de 

4.5 a 10 millimetres avec des yeux rougeatres et un corps verdatre ou vert 
bleute avec des reflets cuivres fait partie de la famille des " Calliphoridae" 
et elle est connue sous le nom "lucilie cuivree australienne" ou tout sim- 
plement " mouches du mouton australien " ou en Latin " Lucilia cuprina" ou 
"Phoenicia cuprina 

En outre, elle a deux paires d’ailes, la premiere paire etant des ailes mem- 
braneuses et la seconde paire etant des ailes posterieures reduites et modifiees 
connues sous le nom " d ’halteres " qui sont utilisees pour la stabilisation du 

vol. 

Le cycle de developpement de cette mouche comprend quatre stades de crois- 
sance : oeuf , larve , pupe et adulte. La lucilie femelle gravide attiree par les 
plaies ou les replis laineux malodorants et humides des moutons pond en 
moyenne 250 oeufs sur la peau de l’animal et qui vont eclore et se muent 
en larves carnivores apres une periode d’incubation ne depasse pas 24//. Ces 
asticots se nourrissent des secretions des plaies et des tissus sous-jacents du 

2 Une myiase designe tout parasitisme d’un etre vivant par les larves de dipteres para¬ 
sites. 
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Fig. 2.2 - La lucilie cuivree. 

mouton pendant trois stades larvaires de duree de 4 a 5 jours. 

Apres la phase larvaire, les larves completement developpees se laissent 
tomber et s’enfoncent dans le sol pour se transformer en pupes donnant des 
nouvelles mouches jeunes. 

Motives par les donnees experimentales obtenues par Nicholson, Gurney , 
Blythe et Nisbet [38] ont propose en 1980, l’equation retardee (2.30) ci-dessous 


qui decrit revolution de la dynamique d’une population au cours du temps 

-SN(t), (2.30) 


dN(t) n / \ / —.V (/.— r) 

w = PN (t - t) exp ' v ' 


dt 


k 


ou 


N (t) represente l’effectif de la population a l’instant t (les lucilies cuivrees 
adultes a l’instant t) 

represente le taux devolution du nombre de la population 
dt 

3 est le maximum de la croissance quotidienne d’oeufs par individu 
k est le nombre maximal d’individus que le milieu peut supported 


3 La capacite de charge du milieu. 
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Mouche femelle gravide attiree 
par les plaies ou les replis 
laineux malodorants et 
humides des moutons 


[ouches adultes 




une odeur caracteristique 
souvent non perdue par 
homme 


Nouvelles 


Environ 250 oeufs pondus 



Mouches 


Les lan es se 
nourrissent de la chair 
du mouton 


Les lan es coinpletement 
developpees se laissent 
tomber au sol 


Les pupes grandissent dans le 
sol pour se transformer en 
nouvelle mouche 


Fig. 2.3 - Cycle de developpement de la lucilie cuivree. 

5 est le taux de mortality par individu (jour -1 ) 

r est la duree de la phase de maturation (le cycle de developpement). 
Plusieurs generalisations et une recherche tres active s’est amorcee recemment 
sur cette equation, par exemple en 2004, Liu et Ge [S3], ont etudie l’equation 
different idle de type neutre du premier ordre suivante : 

y [x (t) — cx (t — T (t))] = -a (t) x(t) + b (; t) x{t-T (t)) e-WW-rW). 

Lib 

En conservant la meme denomination; "equation de Nicholson ", de nou- 
veaux modifications du modele ont ete largement considerees dans l’etude 
des croissances d’autres especes. 

2.5.2 Probleme de deux corps en electrodynamique 
classique 

Dans [28] (voir aussi US ES S3]), Driver a etudie un probleme de deux 
corps en electrodynamique classique avec une interaction retardee. II a consi- 
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dere le mouvement de deux particules chargees se deplagant le long de l’axe 
des x et substitue l’expression pour le champ d’une charge mobile, calculee 
a partir du potentiel de Lienard- Wiechert, dans la loi de la force de Lorentz- 
Abraham. La reaction de rayonnement est omise, mais le retard est inclus en 
raison de la vitesse limitee de la propagation, c, des effets electriques. la mo- 
delisation mathematique de ce probleme conduit a un systeme d’equations 
different idles a retard. En partant de ce modele, il a obtenu un systeme de 
six equations different idles decrivant revolution des etats, des vitesses et des 
retards temporels. 

Pour decrire son modele, nous designons par x t (t) (i — 1,2) les positions 
des deux charges ponctuelles sur l’axe dans un systeme inertiel donne a l’ins- 
tant t. Soit Vi (t) = x\ (t) ( i = 1,2), les vitesses des charges. Comme deja 
mentionne, nous omettons la reaction de rayonnement mais nous permettons 
un champ electrique externe, E ext (t, x) , dans la direction x, qui est suppose 
continu sur un ensemble ouvert I)\ dans le plan (t. x). Alors, l’equation du 
mouvement d’une particule chargee i est donnee par : 

7- , l3 /2 = ^ iE 0 (t’ x i W) + diEext ( t,Xi (t)) , i, j e {1,2} , jV h 

[i - vf (t) /c 2 ] 

(2.31) 

ou m, est la masse au repos et q 1: est la grandeur de la charge i, c est la 
vitesse de la lumiere dans le vide, et E 3 (t, x) est le champ electrique a (t, x) 
du a une autre charge ] ^ i. Le champ magnetique de la charge j n’est pas 
implique dans ce cas unidimensionnel. 

Le champ a l’instant t et au point x produit, par la charge j est suppose 
calcule a partir des potentiels de Lienard- Wiechert. L’expression de ce champ 
implique un decalage temporel t — Tji, representant l’instant ou un signal 
lumineux devrait quitter la charge j pour arriver a x l (t) a l’instant t. Done, 
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le retard doit etre une solution de l’equation fonctionnelle 


Tji (t) = | Xi (t ) - Xj (;t - Tji {t))\ /c. 


(2.32) 


II est clair que Tji (t) ne peut pas etre ecrit explicitement. 


L’introduction des retards dans l’equation (2.31), necessite de determiner 


des trajectoires initiates des deux charges sur un certain intervalle approprie 

[a, t 0 ]. 

Considerons maintenant ces trajectoires initiates et leurs extensions {x\ (t ), X‘> (t)) 
definies sur un certain intervalle [cq/3), ou f3 > t 0 , tel que 

(а) chaque x\ (t) est continue et \x[ (t)| < c pour tout t £ [a, (3) ; 

(б) X 2 ( t ) > X 2 ( t) et ( t , Xi ( t )) £ D pour tout t £ [ t 0 , (3 ); 

(c) les deux equations fonctionnelles r ^ = | x % (t 0 ) — Xj (to — r< ji) | /c ad- 
met des solutions r ( 3, i ^ j g 2} . 


Puis Driver prouve que (xi (t), #2 (t)) est une solution de (2.31) —(2.32) 
si et seulement si elle satisfait au systeme de six equations differentielles a 
retard suivant pour t £ (t 0 ,/3) : 


X[ (t) = Vi (t ), 

(-1 ) l Vi (t) - (-1 Yvjit-Tji (t)) 


T 'n (!) = 


c- i-iyvjit-Tjiit)) 

v'i (t) _ (-1)* CLiC C + (-1)* Vj (■t - Tji (t)) 


, [1 — v i (t) /c 2 ] 


3/2 


T l ($ 


-iy V j(t-Tji (t)) 


+ Qi E ext ( t,Xi (t)) /mi, 
(2.33) 


ou Tji (to) = t^, di = q\qi/ (in eo m, c 3 ) est une constante, et en particulier, 
eo est la constante dielectrique de l’espace libre et (■ i,j ) = (1, 2) ou (2,1). 

Driver [2E] a montre que si les trajectoires initiates donnees satisfont a 
la condition (a) pour a < t < t 0 , la condition ( b ) en t 0 , et la condition (c), 
et si E ext (t,x) est lipschitzienne par rapport a x dans chaque sous-ensemble 
compact de D\ et si la vitesse initiale de chaque particule est lipschitzienne, 
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alors il existe une solution unique. Cette solution pent etre continuee d’exister 
tant que les charges n’entrent pas en collision (lim x\ (t) = lini :c 2 ( t ) quand t 
se rapproche de l’extremite droite de l’intervalle maximal d’existence) et ni 
(t,x i (f)) ni (t,x 2 (f)) ne s’approche de la frontiere de D\. II faut signaler ici 
que dans Particle de Driver et Norris ra. la condition de Lipschitz ci-dessus 
pour les vitesses initiales est relachee a l’integrabilite de la vitesse initiale sur 
[cr, t 0 ]. 

Dans Driver [23], un cas particular a ete donne ou les positions et les vi¬ 
tesses des particules a un instant donne determineront l’etat du systeme. 
Plus precisement, dans cet exemple d’equations electrodynamiques du mou- 
vement, les valeurs instantanees des positions et des vitesses des particules 
determineront leurs trajectoires, si les solutions sont definies pour tous les 
temps posterieurs. Cette propriety etait frequemment conjecturee, affirmee, 
ou implicitement supposee, comme dans la mecanique newtonienne et comme 
indique par la longue liste de references citees dans Particle de Driver [23], 
mais cette propriete ne doit pas etre attendu pour les equations electrody- 
namiques generales. Dans le cas ou E ext (t, x) = 0 pour tout (t, x) G M 2 et 
Q 1 Q 2 > 0 (les deux charges sont de meme signe), lim \x 2 ( t ) — x\ (t)| = 00 et 

t —>-oo 

|Vi (t)| < c pour tout t > a. C’est un resultat tres interessant qui prouve que 
le retard Tji (t) peut devenir non borne, tel qu’on obtient un systeme d’equa¬ 
tions different idles fonctionnelles avec des retards non bornes dependants de 
l’etat. 

Dans le cas ou les mouvements sont consideres tridimensionnels, alors on 
obtient un systeme d’equations differentielles fonctionnelles de type neutre 
avec retards dependants de l’etats, et le taux de changement de u* au temps 
actuel depend aussi de sa valeur historique u' (t — Tji). Plus precisement, si 
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nous introduisons un vecteur unitaire 


Ui = 


Xi — Xj (t — 
GTji 



et une quant ite scalaire 


lij = 1 - ~ V 3 ~ r 3i) ~ U h 

C 


comme dans |23, le point designe le produit scalaire dans M 3 (sans oublier 
que dans ce cas, x \, x 2 sont des vecteurs dans R 3 ), alors la loi de force de 
Lorentz donne 


( t ) = 


qi (l - \vi \ 2 /c 2 ) 


1/2 


rrii 


[Ej + ( Vi/c.Ej ) (iti - Vi/c) - ( Vi/c.Ui ) Ej 


(2.34) 

ou Ej est le champ electrique retarde (valeur vectorielle) arrivant a x, a 
l’instant t de la particule j. Ce champ, dans M 3 , peut etre trouve a partir des 
potentiels de Lienard- Weichert tels que 


E 3 = Wi Tji) /c] [l-\Vj\ 2 (t- Tji )] 

T ji lij 

4 — -^-Ui x ([Ui- Vj (■ t - Tji ) /c] x v'j ('t - Tji)) , (2.35) 

T ji lij 

on k > 0 est une constante en fonction des unites, et x designe le produit 
vectoriel croise dans M 3 . L’adaptation dynamique de t :} l est donnee par 



u,j. [Vi - Vj (t - Tji)] 
c lij 


(2.36) 


D’apres ce qui precede, le mouvement de chaque particule est influence par 
les champs electromagnetiques des autres, et en raison de la vitesse hnie de la 
propagation de ces champs, les equations du modele decrivant le mouvement 
des particules chargees via une action a une distance impliquera des retards 
temporels qui dependent de l’etat du systeme entier. 
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Driver et Hoag [2jB, 33] ont note que si Ton considere que les lois de base 
de la physique sont symetriques par rapport au renversement du temps, done 
l’existence de ces retards implique l’existence des termes avancees dans les 
equations, et done on est conduit a un systeme d’equations different idles 
fonctionnelles avec des arguments mixtes on de type neutre [26] . 

L’heritage des travaux de Driver et ses collaborateurs en electrodyna- 
mique est extraordinaire. Malgre qu’ils ont ete publies depuis 55 ans, ils 
marquent la naissance de la theorie des equations different idles a retard de¬ 
pendant de l’etat et ils ont ete couronnes de nombreux succes dans Inspira¬ 
tion pour le developpement theorique et pratique pour de nouveaux resultats. 

Cependant en pratique, beaucoup d’EDR peuvent deduites a partir d’equa¬ 
tions aux derivees partielles (EDP) ou d’equations aux derivees partielles a 
retard (EDPR). On illustre cela avec trois exemples : revolution au cours 
du temps des populations de cellules souches hematopoietiques, une ligne de 
transmission sans pertes et un modele de reaction-diffusion. 

2.5.3 Hematopoiese 

R ne s’agit pas ici de faire une presentation biologique detaillee de l’he- 
matopoiese, mais plutot de donner une idee de ce qu’est 1’hematopoiese, 
de presenter un modele mathematique pionnier et de citer quelques travaux 
dedies a cette thematique et qui nous paraissent incontournables. L’hemato- 
poiese est un mot venant du grecque haimatopoiein, d’apres poiein "faire" 
et plus precisement, e’est la production et le renouvellement quotidien d’une 
quantite tres importante de cellules sanguines (les globules rouges lors de 
l’erythropoiese, les globules blancs lors de la leucopoiese et les plaquettes 
lors de la thrombocytopoiese) dans les organes hematopoietiques (moelle 
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osseuse chez l’adulte, foie et rate chez l’embryon). 

Cette production des elements figures du sang est assuree par les cel¬ 
lules souches hematopoietiques (CSH) dites primitives et qui ont, la capacite 
d’auto-renouvellement et de differenciation vers l’une ou l’autre des lignees 
cellulaires. Elies se differencient en quatre sous-populations selon leurs fonc- 
tions et leur etat de maturation : cellules souches (immatures), progeniteurs, 
precurseurs (cellules en cours de maturation,) et cellules matures (differen- 
ciees). 


Auto-renouvellement 



Cycle cellulaire 

Taux de proliferation 
Duree de proliferation 

G1/S/G2/M GO 


T 


Apoptose 


Auto-renouvellement 

Differenciation 

Apoptose 


Apoptose 


Maturation 


Mort V 


A, B : Compartments de proliferation 
C, D : Compartments des cellules 
differ enciees 


Fig. 2.4 - Schema representant les differentes etapes de l’hematopoiese [3] 


Les premiers modeles mathematiques remontent a la fin des annees 70 et 
parmi les travaux pionniers qui ferment les bases de l’etude mathematique 
de 1’hematopoiese, on trouve les etudes de M.C. Mackey en 1978 [35] sur 
l’hematopoiese periodique et l’anemie aplasiqucQ 

4 Anemie due a un dysfonctionnement du tissu medullaire, incapable de produire cor- 
rectement les elements figures du sang. 
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Nous allons maintenant presenter le premier modele de Mackey [59] qui 
est inspire par des travaux anterieurs de Lajtha en 1959, Burns et Tannock 
en 1970 et qui decrit revolution au cours du temps de deux populations de 
cellules souches hematopoietiques immatures, une proliferante (les cellules 
sont en phases G\. S, Gb, M) et l’autre quiescente (les cellules sont en phase 
de repos Go). En partant d’un systeme couple d’EDP hyperboliques du pre¬ 
mier ordre et grace a la methode des caracteristiques, il trouva un systeme 
de deux equations differentielles a retard. 



Fig. 2.5 - Le cycle cellulaire [2j. 


On a les equations de transport suivantes : 


d d 

—n (t, a) + —n (t, a) = - (6 + (N (t))) n (t, a) 
qt da 

0 . . o . . , 

di p a ' + da P a ' = ~ 7P a ' 

n ( t , 0) = 2 p (; t, t ) 

P (t, 0) = / 0 + °° P (N (t)) n (■ t , a) da 


a > 0, t > 0 
0<a<r, t > 0 
t > 0 
t > 0, 
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avec les conditions initiates : 

n (0, a) 
P (0, cl) 
et la condition naturelle : 


= n 0 (a) G L 1 ((0, +oo), M + ) 
= Po (a) e L 1 ((0,r) ,M + ), 


(2.38) 


ou 


a : 


n ( t , a) 

P (t, a) 

N (t) : 

5 : 

7 : 

P- 

T : 


lim n (t, a) = 0, t > 0, 

-t+CO 


Page qui est le temps passe par une cellule dans une des 
deux phases, quiescence ou proliferation, 
la densite des CSH en phase de repos a l’instant 
t > 0 et d’age a > 0. 

la densite des CSH en phase de proliferation a 
l’instant t > 0 et d’age 0 < a < r. 

la densite de la population totale des cellules dans la phase de repos 
ou N (t) = / 0 + °° n (t, a) da. 

Le taux de mortality des cellules au repos (jour -1 ), 
le taux d’apoptose (mort naturelle, suicide cellulaire, ou mort 
cellulaire programmee) (jour -1 ), 
le taux de reintroduction (jour -1 ). 

la duree moyenne de la phase de proliferation, (la duree du 
cycle cellulaire) (jour). 


Dans le systeme (2.37), la premiere condition frontiere signifie que toutes les 


cellules proliferantes entrent immediatement apres leur division dans la phase 
Go. La deuxieme condition frontiere decrit le flux des cellules quiescentes qui 
entrent dans la phase de proliferation. 

En vertu de la methode des caracteristiques, on arrive a 

e -7 V (0, r — t) 0 < t < t 


P (t, r) = 


e 1T p (t — r, 0) t > 0. 


(2.39) 
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En posant 


P (t) = / p (t, a) da, 


(2.40) 


la densite des cellules dans la phase de proliferation a l’instant t > 0 et en 


tenant compte des conditions frontieres et de (2.39), l’integration du sys¬ 


teme (2.37) par rapport a Page donne le systeme couple de deux equations 
differentielles a retard temporel discret suivant : 


I N'(t) = -(5 + f3(N(t)))N(t) + 2e~^ T f3(N(t-T))N(t-T) t>r 

\ p> (t) = -jP (t) + (3 (N (t)) N (t) - e~ 1T /3 (N (t - r)) N (t - r) t> r, 

(2.41) 

avec des conditions initiates pour t G [0, r]. 

Dans la premiere equation, la perte de la non proliferation des cellules a la 
phase de proliferation et a la differentiation est representee par le premier 
terme de droite tandis que le deuxieme terme represente la phase de produc¬ 
tion des cellules dans la phase de repos G 0 a partir de la proliferation des 
cellules souches. Le facteur 2 compte pour Famplification des effets de la di¬ 
vision cellulaire tandis que e~ 1T compte pour F attenuation dans la phase de 
proliferation en raison de Fapoptose. Ce modele a ete modifie de nombreuses 
fois. 


2.5.4 Ligne de transmission sans pertes 

Maintenant, en vertu de la transformation de d’Alembert, on va voir com¬ 
ment pent on reecrire une EDP hyperbolique en tant que systeme a retard 
de type neutre m similaire a [49]. 

Cet exemple traite une ligne de transmission sans pertes, a une extremite 
(x = 0), il y a une source de tension constante E et a Fautre extremite 
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(x = /), une capacite est branchee en parallele avec une diode tunnej^} 


R x=0 x=f 



1 


id 


Fig. 2.6 - Circuit de la ligne de transmission ra 

Le courant i(.,.) et la tension v (.,.) sont des fonctions de t et x qui satisfont 
l’equation du telegraphe, done on obtient les deux equations aux derivees 
partielles du premier ordre de type hyperbolique suivantes : 

r di dv „ ^dv di „ , . . 

L~dt dr = ^ ^ ~dt ~dr = ^ ® — x — ^ (2.42) 

avec les conditions aux limites 

dv 

E = Ri(t, 0) + u(£,0) et i(t,l) = Q — (t,l) + i d (v(t, /)), (2.43) 

ou l est la longueur de la ligne, L et C sont l’inductance et la capacite du 
conducteur par unite de longueur, R est la resistance a l’entree, C'i est la 
capacite en parallele avec la diode tunnel, id(v) est la courbe caracteristique 
courant-tension de la diode (une fonction polynomiale non lineaire) donnee. 
Les points d’equilibre (ty, to) satisfont une equation lineaire donnee par 

E = v 0 + Rio, io = id {vq) ■ 

5 Ce genre de diodes est utilise dans les amplificateurs haute frequence de circuits os- 
cillants electroniques par exemple. 
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-0.0004 


0.0016 
0.0014 
0.001 


0.0008 
0 0006 
0.0004 
0.0002 


Fig. 2.7 - Courbe caracteristique courant-tension et ligne de charge de la 

diode tunnel. 133 


On suppose que le point de fonctionnement est donne par 1’approximation 
du premier ordre 

i(v d ) « i(v 0 ) + m (v d - v 0 ), 


ou m est une constante positive. La formule de D’Alembert et les conditions 


aux limites (2.43) nous permettent d’obtenir l’equation differentielle lineaire 
de type neutre suivante : 


d^> 

~di 




« - T) = A<t(Z) + B'I'R 


T ), 


avec A 



r = 2\jLCl et ^ =: t — bl. 



(2.44) 
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Ce qui prouve une transformation du systeme de deux EDP de type hyper- 


bolique (2.43) en une equation a retard de type neutre. 


2.5.5 Equation de reaction-diffusion avec argument de- 


vie 


L. Sergeeva [76] . A.M. Samoilink et collab [72] [73] ont etabli les conditions 
dans lesquelles il est possible de construire des solutions globales d’une EDP 
avec argument devie dans la variable temporel. Ils ont montre l’existence et 
l’unicite de ces solutions en vertu du theoreme de point fixe de Banach et 
ils ont donne la structure de ces solutions et l’algorithme utilisee pour leur 
construction par Futilisation de la methode de separation des variables et la 
methode des suites recurrentes. 

Considerons le probleme aux derivees partielles lineaire non homogene 
avec argument devie suivant : 

u t (x, t) = p (t) u xx (x, t + p) + q (x, t ), (x, t) e Q, 

u (0, t) = u (l, t) = 0, t e M, 

avec des conditions aux limites de Dirichlet homogenes 


(2.45) 


u (0, t) = u (/, t) = 0, t G M, 
ou Q = {(x, t) : 0 < x < /, t e M} et p e M\ {0} . 


Remarque 2.2 Dans le cas on t > 0 (la solution n’existe en general que 


dans le futur), l’EDP avec argument devie (2.45) est de type parabolique et 


dite " equation de reaction-diffusion avec argument devie " ou " equation de 
la chaleur non homogene avec argument devie " |otm (x,t) est une fonction 
G La transformation n’est pas unique. 

' Le premier exemple a ete etudie par Fourier en 1822 dans son oeuvre revolutionnaire 
" Theorie analytique de la chaleur". 
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definie pour la position x et le temps t et represente les individus qui peuvent 
etre differents d’un probleme a l’autre; ils peuvent representer des concen¬ 
trations, des densites, des charges electriques, des temperatures, etc. La pre¬ 
miere derivee par rapport au temps traduit la vitesse de formation, q(x,t) 
est une fonction appelee "la reaction", "le terme reactif" ou "la source " et 
elle decoule de toute interaction entre les individus, p (t) est le coefficient de 
diffusion et p ( t ) u xx ( x , t + p) est le terme de diffusion qui peut representer la 
diffusion moleculaire ou un certain mouvement aleatoire d’individus dans une 
population et /i est l’argument devie. Les conditions aux limites homogenes 
de Dirichlet signifient qu’il n’y a pas d’individus sur la frontiere. 


Pour simplifier l’EDP avec argument devie (2.45), Samoilenko et all m 
CUES], ont utilise la methode de separation des variables sur le probleme 
homogene qui a comme fonctions propres, 

knx 


X k (x) = sin 


l ’ 


associees aux valeurs propres 


7 rk 


Afc = [ — ] , k — 1,..., n pour certain n > 1. 


II ont construit une solution globale u (x, t ) du probleme (|2. 45) sous la forme 

(2.46) 


knx 


u[x,t) = y^T k (t) sin ——, ( x,t)eQ. 


k= i 


l 


L. Sergeeva |76j a suppose aussi que le terme source q (x. t ) peut etre re¬ 
presente comme une somme des n premiers termes de la serie de Fourier 
impaire 


q (x, t ) = ^qk (t) sin 


k =1 


knx 

T’ 


ou 


Qk {t) = J [q k (^, t) sin 


l 


(2.47) 


(2.48) 
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Done, la substitution de (2.47) et (2.48) dans (2.45) donne l’equation a ar¬ 
gument devie suivante : 


T' k ( t ) + X k p ( t ) T k (t + n) - q k ( t ) = 0, k = 1,..., n. (2.49) 


2.6 Recapitulatif 

Revenons un instant sur les principaux points abordes dans ce chapitre. 
Premierement, on a presente un bref historique sur les equations fonction- 
nelles a retard suivi par une classification des differents types de ces equations 
cites dans la litterature. Ensuite, on a propose un rapide survol sur la theo- 
rie d’existence et d’unicite des EDR et des EDN et quelques methodes de 
resolution comme la methode des etapes, celle de Runge-Kuta et la transfor¬ 
mation de Laplace. Finalement, on a expose cinq modeles d’equations avec 
retard, les deux premiers modeles (le modele de Nicholson et le modele de 
deux corps en electrodynamique classique) peuvent etre modelises mathe- 
matiquement par des equations differentielles a retard discret ou variable ou 
bien par equations differentielles de type neutre tandis que le troisieme et 
le quatrieme phenomenes sont modelises par des EDP sans retard mais leur 
etude mathematique conduit a obtenir des EDR ou des EDN. En effet, dans 
le modele decrivant Involution au cours du temps des populations de cellules 
souches hematopoietiques, la methode des caracteristiques conduit a obtenir 
un systeme couple de deux equations differentielles a retard temporel discret 
et dans le modele d’une ligne de transmission sans pertes, la transformation 
de d’Alembert nous permet de reecrire l’EDP hyperbolique en tant qu’une 
equation differentielle a retard de type neutre. 

Enfin, on a presente un modele d’une equation aux derivees partielles 
lineaire non homogene a argument devie ou la methode de separation des 
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variables sur le probleme homogene conduit a une equation differentielle a 
argument devie. 

Pour plus de details nous referons aux livres de Hale [3T|, Hale et Lunel 
[42J et le lecteur pourra consulter aussi le livre de Kolmanovskii et Nosov 
ra pour des exemples de modeles comportant des termes a retard. 
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.Existence et positivite des solutions periodiques d’une 
classe d’equations differentielles fonctionnelles de type 

neutre du second ordre 


Sommaire 


13.1 Introduction! 

13.2 Preliminairesl 


3.3 Existence de solutions periodiques positives 

3.4 Recapitulate . 


62 

63 

72 

78 


C e chapitre est consacre a prouver l’existence de solutions periodiques 
positives d’une classe d’equations differentielles du second ordre non 
lineaires de type neutre. L’etude est faite en vertu d’une approche basee sur 
1’application du theoreme de Krasnoselskii combine avec quelques proprietes 
de la fonction de Green d’une EDO lineaire du second ordre. 
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Chapitre 3. Existence et positivite des solutions periodiques d’une classe 
d’equations difFerentielles fonctionnelles de type neutre du second ordre 

3.1 Introduction 

Dans ce chapitre, nous sommes essentiellement inter esses par l’etude de 
l’existence de solutions periodiques positives d’une classe d’equations diffe¬ 
rent ielles non lineaires du second ordre de type neutre qui apparaissent dans 
la modelisation de nombreux phenomenes rencontres en physique, biologie, 
mecanique, economic, etc. 

On peut considerer que l’etude des equations fonctionnelles de type neutre 
commence avec les trois article^ (presque) oublies de I.P. Kabakov (1946), 
Sokov (1946), I.P. Kabakov et A.A. Sokov (1946) sur les phenomenes de pro¬ 
pagation. A partir des annees 60, ces equations ont ete l’objet de plusieurs 
travaux. Notamment, on trouve les travaux de A. Bolinia et Myshkis (1960), 
Brayton et Mille -Luxriante (1964), Cooke et Krumme (1968), Cooke (1970), 
Halanay et Rasvan (1970) et Filimonov (1981), etc. 

Durant ces deux dernieres decennies, de nombreux travaux ont etabli plu¬ 
sieurs resultats sur l’etude des solutions periodiques. Certains mathemati- 
ciens ont utilise des transformations pour se ramener a une equation plus 
simple on a un systeme d’equations. D’autres chercheurs ont donne la solu¬ 
tion sous forme de serie convergente vers la solution exacte et certains d’entre 
eux ont traite les equations fonctionnelles du second ordre de type neutre par 
des methodes numeriques telles que la methode de Ritz , la methode des diffe¬ 
rences finies, la methode des elements finis, la methode des splines cubiques 
et la methode des multi-derivees. 

Dans ce chapitre, ces methodes habituelles peuvent sembler inefficaces pour 
traiter l’equation differentielle non lineaire du second ordre de type neutre 

1 Les trois articles ont ete ecrits probablement plusieurs annees plus tot rnais publies en 
1946 . 


62 





Chapitre 3. Existence et positivite des solutions periodiques d’une classe 
d’equations difFerentielles fonctionnelles de type neutre du second ordre 


suivante : 

^x(t)+p(t)^x(t) + q(t)x(t) = ^Q(t,x(t-r(t))) (3.1) 

+f(t, hi (. x(t )), gi (x(t - r(t)))), 

ou p et q sont des fonctions positives continues. La fonction Q : IxR —» M 
est differentiable, hi, gi : K —> Eet/:lxlxK —* M sont des fonctions 
continues par rapport a leurs arguments. 

On va utiliser le theoreme de point fixe de Krasnoselskii combine avec quelques 
proprietes d’une fonctions de Green pour etablir les resultats desires. L’idee 
de l’approche est de convertir l’equation consideree en une equation integrate 
qui s’adapte a l’application du theoreme de point fixe de Krasnoselskii pour 
prouver l’existence de solutions periodiques positives. 

3.2 Preliminaires 

Dans ce paragraphe on donne quelques definitions et certains resultats essen- 
tiels qui representent un outil important pour le reste de ce chapitre. 

Soit 

Pt = {x (t) G C (M, M) : x (t + T) — x (t), T > 0} , 
muni de la norme 


||x||oo = sup |x(f)| = sup |x(i)|. 
teR te[o,T] 

Alors (Pt, ||-||oo) es t un espace de Banach. 

Puisque nous sommes interesses par l’existence des solutions periodiques de 


l’equation (3.1), il est naturel de supposer que 


p(t + T) = p(t), Q(t + T) = q(t ) et r(t + T) = r(t), (3.2) 
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ou r est une fonction continue, et r(t) ^ r* > 0. En outre, nous supposons 
que 

/ p(s)ds >0, / q(s)ds > 0. (3.3) 

J 0 Jo 

Nous supposons aussi que les fonctions Q(t, x) et fit, x. y ) sont T—periodiques 
en t. ce qui est equivalent a 

Q[t + T,x) = Q(t , x) et f(t + T, x, y ) = f(t, x, y). (3.4) 

D’abord, on considere l’equation suivante : 

x' (■ t ) + a (t) x ( t ) = AM (x (t)) + m{t) , (3.5) 

ou a est une fonction continue T— periodique et non negative avec 

/ a(s)ds > 0, 

Jo 

m est une fonction continue T—periodique, A > 0 est un parametre et 
M : M + —i M + est une fonction continue. 


Lemme 3.1 J55 1/ Supposons que 


— N 6 (0, +co]. 


(3.6) 


x —>+oo X 

Alors, {3.5 ) admet au moins une solution T—periodique positive si £i < £ 
et A G (£i,£ 2 ] , ou ^ et sont definies comme suit : 


£i = 


exp ( J q T a(u)duj — 1 


Rn 


NaT 


') S 2 — 


exp ( J q T a(u)duj — 1 


QiT 


avec 


Ro = \\R 


Q i = max V (x ), 

x&[0,Ro(a+l/<j)\ 

r t+T 


a = exp ( — / Li(u)du ) , 


V{x) = 


R(t) — / G\ [t, s ) m (s) ds > 0, G\ (t, s ) = 


M (x) si x G [0, +oo) 

M (0) si x E (—oo, 0) 

exp a{u)du) 


exp 


J q T a{u)du)j 
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Lemme 3.2 /55j/ Supposons que (3.2) et (3.3) sont verifiees et 


Ri 


exp ( Jq p(u)du ) — 1 


QiT 


> 1 , 


(3.7) 


ou 


et 


R\ = max 

te[o,T] 


r t+T exp ( I p(u)du 

r -— q(s)ds 

exp ( / p(u)du j — 1 


Qi=^l + exp^y p^di^j^j R\. 

Alors, il existe des fonctions continues T—periodiques a et b telles que 
b(t) > 0, Jq a(u)du > 0 et 

a(t ) + b(t) = pit), --b(t ) + a(t)b(t ) = q(t) pour tout i G R. (3.8) 

C lb 


Preuve. De (3.8), on obtient 


b' ( t)+p(t)b(t ) = b 2 (t) + q(t). 


(3.9) 


Soit Mix) = x 2 , m(t) = q(t) et A = 1. Alors 

AT 1- Mix) 

N = lim -- = +oo. 

x —>+oo x 


Par le Lemme 3.1, si 


R 0 [exp( / p(u)du) — 1] 

0 < A = 1< -—- 

QiT 


alors, (3.9) admet au moms une solution positive T— periodique continue bit). 
Maintenant, comme 

/ q{s)ds > 0, 
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et 


il vient 


b ' {t) ' a(t)= ,W 


b(t) 


Kty 


a(u)du = [ \du > 0. 

) Jo Ku) 


D’ou le resultat desire. 


Lemme 3.3 [8U] Supposons que les conditions du Lemme 3.2 sont verifiees 
et $ G Pt, alors Vequation 


d 2 


d 


jpx{t)+p{t)—x{t) + q{t)x{t) = $(t), 


(3.10) 


admet une solution T—periodique. En outre, toute solution de (3.10) s’ecrit 
sous la forme : 

r t+T 

(3.11) 


x(t) — J G(t, s)Q(s)ds, 


ou 


G(t,s ) 


f( exp [ [(' b{v)dv + a{y)dv\ du + f* +T exp 

f t u b(v)dv + f* +T a(v)dv 

du 


exp ^ a{u)duj — 1 

exp ^ b(u)duj — 1 



(3.12) 


Preuve. Definissons E a et Eh comme suit : 

E a = exp I / a (u) du ) — 1, Eh = exp I / b (u) du ) — 1. 


Par un calcul direct, nous pouvons voir que (3.11) est une solution T—periodique 


de (3.10). 


Maintenant, on suppose que x(t) est une solution T —periodique de (3.10), 

nous avons 


en vertu du Lemme 3.2 


^z(f) + + a (t)b(t)x(t) = $(t), 
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ce qui conduit a 


d 


—x(t) + b(t)x(t) = 


rt+rp exp I / a(u)du 


exp ( I a(u)du ) — 1 


-<E>(s)ds. 


Done 

t+T exp I / b(u)du 
x(t) = — U 


rs+T exp I / a(u)du 


-$(u)du 


E a E b J t 


exp yj a(u)du J — 1 

rs+T / r v \ 

exp I / a(u)du] $(u)duds 


E a E b J t 


E a E b J t+T Jv-t 


exp yj b(u)duj — 1 

rt+T / rs 

/ exp I / b(u)du 

t dv^ exp b(u)dv^j exp a(u)du^j $(u)ds 

t+2T rt+T / rs \ f f V \ 

dv j exp I / b(u)du ) exp I / a(u)du I $(u)ds 

1 rt+T rs / ru rs \ 

$(s)ds / exp I / b(v)dv + / a(u)du J du 

t+T rt+T / ru rs+T \ 

&(s)ds / exp I / b(v)dv + / a(v)dv ) du 

Jt \Jt Ju J 


ds 


E a E b J t 


E a E b 

rt+T 


G(t , s)<l?(s)ds. 

J t 

Ce qui etablit le resultat desire. 


Corollaire 3.1 [8W La fonction de Green G deftnie par (3.12) satisfait les 
proprietes suivantes : 

G(t, t + T) = G(t, t ), G{t + T, s + T) = G(t , s), 

o exp( / 6(u) du) 

—G(t,s) = a(s)G(t, s) - JjJ -, 


exp( / b ( v )du) — 1 

Jo s 

exp( / a (u) du) 


—G(t,s) — —b(t)G(t, s) H- -4jJ- 


exp( / a (u) du) — 1 

Jo 
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Remarque 3.1 Par un calcul simple, on obtient que la derivee seconde par 
rapport a s de la fonction de Green G(t. s ) est comme suit : 

—G(t, s) = (a 2 (s) + a'(s) j G(t, s) - p (s) — 

ex 

Le Lemme suivant est important pour etablir l’existence. 


exp( / b (■ v ) dv) 


p( / b (v) dv) — 1 


Lemme 3.4 Supposons que les conditions (3.2)-(3.4) et (3.7) sont satis- 


faites. Si x € Pt, alors x est une solution de I’equation (3.1) si et seulement 
si 

pt +T 

x(t) = Q(t,x(t-T(t)))- p(s)E(t,s)Q(s,x((s~r(s)))ds 

nt +T 


+ 


G(t, s ) [/ (s, ht (x (s)), g\ (x (s — t (s)))) 


(3.13) 


+ ( a'(s) + a 2 (s )) Q (s, x((s - r (s))) 


ds, 


ou 


Preuve. 


E(t,s) 


exp(// b (y) dv) 
exp(J Q T b (v) dv) — 1 


(3.14) 


Ligne directrice : Par Vutilisation du Lemme 3.3 et certaines proprietes 
de la fonction de Green citees dans le Corollaire 3.1, on peut convertir 
I’equation (3.1) en I’equation integrate (3.13). 


Soit x G Pt une solution de l’equation (3.1). D’apres le Lemme 3.3, on 


a 


x(t) = J G(t, s) ^Q(s,x((s-t(s))) (3.15) 

+/ (s, hi (x (s)), g! (x (s — r (s))))] ds. 
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Une integration par parties nous donne 


rt+T 


G(t, X (( S ~ T ( S ))) ds 


G(t,s)-^Q(s,x((s-T(s))) 


t+T 


ft +T Q Q 

s) —Q (s, x((s - t (s))) ds. 


Comme, 


G(t, s)-^Q {s, x ((s - t (s))) 


-i t+T 


= 0 , 


alors, 

r t +T 


d 2 


rt +T 


d 


d 


G(t, s )g^Q (s, x((s - t (s))) ds = - I -G(t,s)—Q(s,x((s-T(s)))ds. 


ds 


ds 


Par une deuxieme integration par parties, on obtient 


ft+T 


G(t, s )j^Q (s, x((s - r (s))) ds 
Q(s,x((s-T(s)))^G(t,s) 


ft+T 


d 2 


+ Q(s,x((s-T(s)))—G(t,s)ds. 


Posons 


I = 


Q(s,x((s-T(s)))-^G(t,s) 


t+T 


J t 


Par le Corollaire 3.1 


I = 


( 


Q(s,x((s - r(s))) 


a(s)G(t, s ) 


exp( / b (v) dv) 


\ 


-i t+T 


\ 


exp (f b (v) dv) — 1 


/ J 
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done 


I = -Q(t + T,x({t + T-r{t + T ))) 


/ 


ct+T 


X 


exp( / b (v) dv ) 


a(t + T)G(t , t + T) 


V 


+Q (t, x((t-T (t))) 


ex 


P(/o b ( v ) dv ) - 1 


/ 


( f* \ 

exp( / 6 (u) du) 

a(t)G(t,t )- / pT * - 


V 


ex 


P(/o & (u) *0 - 1 


i.e., 


/ 


/ = -Q(/..r((/ -r(/))) 


exp( / 6 (n) du) 




a(t)G(t,t ) — 


V 

+Q (t, m((t - r (£))) f 'a(t)G(t, t) 


exp (f b (v) dv) — 1 
1 


exp (f b (v) dv) — 1 


Comine 

ft +T 


(&_ 


j \ J Q( s ^ x (( s ~ r ( s ))) ds 

ft +T 

— j { (a(s) + a 2 (s )) Q (s, x((s - r (s))) G(t, s) - p (s) E (t, s) Q (s, x((s 
Alors, 

ft +T / Q2 \ 

J y-g^G(t,s)jQ(s,x((s-T(s)))ds 
= Q(t,x((t- r(t))) 

pt +T 

+ / {(a(s) + a 2 (s)) Q (s,m((s - r (s))) G(t, s) - p (s) E ( t,s) Q (s,x((s - 

(3.16) 


ou E est donnee par (3.14). Enfin la substitution de (3.16) dans (3.15), achieve 
la preuve. ■ 


r(s)))}ds. 


r(s)))}ds, 
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Lemme 3.5 {8Bj Soit 


A = p(u ) du, B = T 2 exp 



T 


In (q (u)) du I . 


Si 


A 2 ^ 4B, 


(3.17) 


alors , on a 


(u) du, J ^ — VA 2 — 4 B'j l, 


mm < j a 

max | J a(u)du,J b (u) duX ^ - 


4 B ) := m. 


Preuve. Comme a (t) + b ( t ) = p(t) , b' ( t ) + a ( t ) b ( t ) — q(t), on a 


/*T pT pT pT pT q /^\ 

/ a(u) du+ b(u)du= p (u) duet a (u) du= . du. 

/o Jo Jo Jo Jo h W 


Grace a l’inegalite 

f b 

/ P{ x ) f ( x ) dx 


A exp 


/ f T \ 

/ p (x) In / (x) dx 

I o 


p (x) dx 


i; 


p (x) dx 


ou p (x) ^ 0, / p (x) dx > 0, inf / (x) > 0, nous arrivons a 

J a 

1 f T /If 51 ,, ) 1 f T q(u ) /1 f T q(u) 

— b(u) du ^ exp ( — / In b (u) du I , — / du ^ exp ( — / In ; 

TJo 1 ' P Vr7o ' / T J 0 b (u) " \Tj 0 b(u ) 


Done, nous obtenons aisement 


pT pT pT pT q /^\ 

/ a(u) du b (u) du = b (u) du / du 

Jo Jo Jo Jo b{u) 

^ T 2 exp^^y' In q{u)du 

f T f T f T 

Sachant que / a (u) du + / b(u) du = / p ( u ) du, la conclusion decoule 

Jo Jo Jo 

immediatement de ce qui precede. ■ 


du 


71 














Chapitre 3. Existence et positivite des solutions periodiques d’une classe 
d’equations difFerentielles fonctionnelles de type neutre du second ordre 


Corollaire 3.2 


Les fonctions G et E veriftent 

■T 


i y 


< G (t, s ) ^ 


T exp (^J p(u) dv^j 


(e l - iy 


|£7 (t, s)[ < 


e l — 1 


(3.18) 


3.3 Existence de solutions periodiques posi¬ 
tives 

Pour appliquer le theoreme de point fixe de Krasnoselskii , nous avons besoin 
de definir un espace de Banach (B, ||.||), un sous ensemble convexe ferine © 
de B et de construire deux applications, l’une est une contraction et l’autre 
est continue et compacte. 

Soit (B, ||.|| OQ ) = (Pt, (MU et © = {p e B : K ^ p ^ L}, ou K est une 
constante non negative et L est une constante positive. Nous exprimons 


l’equation (3.13) comme suit 


p (t) = (»¥>) (t) + (Ap) (t) := ( Up ) (t) , 


ou A, B : 


sont deux applications definies par 


rt +T 


(Ap) (t) = j G(t,s)[f(s,hi(p(s)),g 1 (p(s-r(s)))) (3.19) 

+ fa'(s) + a 2 (s)) Q (s, p((s - r (s))) ds, 


et 


(Bp)(t) = Q(t,p(t-r(t ))) 

- J V (- s ) E (t, s) Q (s, p((s — t (s))) ds. 


(3.20) 
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Pour simplifier les notations, on introduit les constantes suivantes : 


T exp( / p ( u ) du) 


a = 


(e* - 1)' 
6 = max \b(t )\, 

te[o,T] 


P = 


e l — 1 


7 = 


(e m - 1) 


2 > 


(3.21) 


u — min |p(t)|, A= max |p(£)| 
te[o,T] ' te[o,T ] 1/1 


On suppose que la fonction Q(t,x ) est localement lipschitzienne et continue 
par rapport a la variable x. done il existe une constante strictement, positive 
k telle que : 


| Q(t,x) - Q(t,y) | < k ||x - y\\, \/t e [0,T] ,Vx e ID. 


(3.22) 


Pour etablir l’existence des solutions periodiques positives de (3.1), on dis 


tingue deux cas Q(t, x) > 0 et Q(t, x) < 0 pour tout t G M et x G B. 
Commenqons d’abord par le premier cas ou Q(t. x) > 0. 

On suppose qu’il existe une constante non negative k\ et deux constantes 
positives /C 2 et a telles que : 


E (t, s) > a, V (t, s ) G [0, T] x [0, T ), 
k\X < Q(t, x) < k 2 x, Vt G [0, T] , Vx G ID, 


et 


k 2 < 1, 

ainsi, pour tout s G [0, T \, x, y G B 

(1 - fci) K + A (3k 2 TL 
7 T 


(3.23) 

(3.24) 


(3.25) 


< f (s.^^x) ,g l (y))+(a(s) + a 2 (s)^jQ(s,y) (3.26) 

(1 — k 2 ) L + jicrkiTK 


< 


aT 


Lemme 3.6 Supposons que les conditions (3.2)-(3.4), (3.7), (3.17) et (3.23) 


(3.26) sont satisfaites, alors I’operateur A defini par (3.19) est compact. 
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Preuve. 


Ligne directrice : En suivant les memes demarches que celles 
utilisees dans I’exemple 1.6 et en vertu du theoreme d’Ascoli — Arzila, 
on peut montrer la compacite de l’operateur A. 


Soil A l’operateur donne par (3.19). II est clair que A est continu, ainsi 


on peut montrer aisement que (A<p) (t + T) = (Aip) (t) 
Pour t G [0, T] et ip, ip E B, on a 

r t +T 


\(A<p)(t)\ < 


G(t, s ) [/ (s, h x (<p (s)) ,g 1 {(p{s-T (s)))) 


ds 


< aT 


+ (a (s) + a 2 (s)J Q (s, <p((s - r (s))) 
(1 — k- 2 ) L + ficrkiTK 


aT 


= (1 — k 2 ) L + fiakiTK. 


D’ou 


Halloo < (1 - k 2 ) L + fiakxTK, 
et par consequent, l’ensemble A (B) est uniformement borne. 

Montrons maintenant que A (B) est equicontinu. Pour cela, soit ip n G B, avec 
n un entier naturel. On va calculer — (A^ n ) (t). puis on va montrer que cette 


derivee est uniformement bornee. A l’aide de (3.2), (3.3) et (3.4), si on prend 


la derivee par rapport a t dans (|3.19[). on obtient 

exp (f* a ( v ) dv) 


d 


S (A>„)M = 


- b(t)G(t , s) + 


exp ^ J o r a (■ v ) dv) — 1 j 
x [/ (s, h x (<p n (s)), g 1 (<p n (s - t (s)))) 

—a (s) Q (s, (p n ((s — t (s)))]. 


En utilisant (3.21) et (3.26), on arrive a 


Jt (i) 


<T(» a + A l ~ k2)L + ,laklTK <D, 

aT 
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pour certaine constante positive D. Alors la suite (Ap n ) est equicontinue. 
Le theoreme d 'Ascoli-Arzela assure l’existence d’une sous suite (A<p n ) de 
(Ap n ) qui converge vers une fonction T—periodique et continue. Par conse¬ 
quent A est continu et A (B) est contenu dans un sous-ensemble compact de 

B. ■ 


Lemme 3.7 Supposons que la condition (3.22) est satisfaite. Si 


(1 + A j3T)k< 1, 


(3.27) 


alors I’operateur B defini par (3.20) est une contraction. 


Preuve. Soit B l’operateur donne par (3.20). II est clair que Bp est continu 
et il est aise de voir que 


(Bp) (t + T) = (Bp) (t ). 

En outre, pour tout t G [0, T] et ip, ip G B, on a 

10M (0-W) (t) I 

< \Q{t,<p(t-T (t))) -Q(t,i>(t-r(t)))\ 

+J p(s)E(t,s)\Q(s,<p(s-T(s))) -Q(s,ip(s-T(s)))\ds 

< (1 + \j3T) k \\p - V'lloo , 

il s’ensuit alors que 


\\Bp - BipW^ < (1 + A/IT) k \\<p - ^|| c 


Nous concluions par la condition (3.27) que B est une contraction. 


L’equation integrate obtenue (3.13) est une somme de deux applications; 


une contraction et une application compacte et continue, done dans ce cas 
on va appliquer le theoreme de point fixe de Krasnoselskii. 
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Theoreme 3.1 Supposons que les conditions (3.2)-(3.S), (3.1), 3.11 et (3. 


(3.21) sont satisfaites, alors Vequation (3.1) admet au moins une solution 


positive T—periodique dans Vensemble B. 

Preuve. 


Ligne directrice : La compacite de A et la contractivite de B sont 


assurees par les Lemmes 3.Q et \3.7\ done pour appliquer le theoreme 
de Krasnoselskii, il nous reste qu’a verifier sa premiere condition. 


D’apres le Lemme 3.6 , A : B —> B est compact et continu et d’apres le 
Lemme |3.7[ B : B —» B est une contraction. De plus, si <p, ^ G B, nous 


pouvons voir que 

(A<p) ( t ) + {Bi>) ( t ) 

= Q (t,ip (t — r (t))) — J p (s) E (t, s) Q {s — t (s))) ds 

+£ G(t, s ) [f (s, h x (p (s)) ,gi(<p{s-r (s)))) 

.2 !„\ \ .-<( „ - (~w\ d s 


+ (a (s) + a (s)J Q (s, <p((s - r (s))) 

nt +T 

< k 2 L — fia / Q (s, ip((s — t (s))) ds 


et +T 


-\-ot 


[f (s,h 1 {p (s)),pi (p (s — T (s)))) 


ds 


+ (a (s) + a (s)) Q (s, p((s - r (s))) 


z , r , (1 - k 2 )L + pak l TK 

< kpL — fiakiTK + a -—-= L, 

aT 


et 


{■Aip) (t) + W) (t) 

f*t +T 

> k\K — X/3 Q (s,p((s — r (s))) ds 

f*t +T 

I [f (s, hi (p (s)) ,gi(ip(s-T (s)))) - a(s)Q (s, p({s - r (s)))] ds, 
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par consequent 

(Aip) (t) + (Bip) (t) > kiK — \f3k 2 TL + 7 
> K. 


(l-k 1 )K + \pk 2 TL 
7 T 


Maintenant, toutes les hypotheses du Theoreme de Krasnoselskii |1.5| sont 
satisfaites, alors il existe un point fixe x G B de l’operateur 


H = A + B, 


i.e., 


x = Ax + Bx. 


Enfin, le Lemme 3.4 assure que ce point fixe est une solution de (3.1). ■ 
Dans le cas ou 0{t. x) < 0, on suppose qu’il existe une constante negative k 3 
et deux constantes non-positives Aq et a telles que : 

k 3 x < Q(t,x) < k±x, Wt G [0, T\ , Vx G ID, (3.28) 

- k 3 \(3T < 1, (3.29) 

ainsi, pour tout s G [0, T\ , x, y G B 

K (1 + Aq/icrT) — k 3 L 
7 T 

< fis,h 1 (x) 1 g 1 (y))+^a(s) + a 2 (s)^Q(s,y) (3.30) 

(1 + hXfiT) L - UK 


< 


aT 


Theoreme 3.2 Supposons que les conditions (3.2)-(3.3), (3.1), 3.11 , (3. 


et ( 3.21)-(3.30 ) sont veriftees, alors Vequation (3.1) admet au moins une 
solution positive T—periodique dans Vensemble B. 


Preuve. La preuve est analogue a la preuve du Theoreme 3.1 
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3.4 Recapitulatif 

Dans ce chapitre, nous avons aborde la question de l’existence des solu¬ 
tions periodiques positives d’une classe d’equations differentielles non li- 
neaires de type neutre et nous avons obtenu des resultats d’existence pouvant 
s’etendre aux EDN d’ordre superieur a deux. Avant d’appliquer le theoreme 
de point fixe de Krasnoselskii pour prouver l’existence des solutions, nous 
avons converti l’equation differentielle de type neutre en une equation in¬ 
tegral s’ecrivant comme somme de deux applications; la premiere est une 
contraction et l’autre est continue et est prouvee compacte moyennant le 
theoreme d 'Ascoli-Arzela. 


78 




CHAPITRE 4 _ 

.Existence et unicite des solutions periodiques d’une 
classe d’equations differentielles iteratives du second 

ordre 


Sommaire 


|4J 

RL2 - 

Introduction!. 

Preliminairesl. 

4.3 

Existence de solutions periodiques. 



4.4 

Existence et unicite des solutions periodiques . 

4.5 

Recapitulate]. 


80 

83 

85 

94 

95 


D ans le but de prouver l’existence de solutions periodiques d’une classe 
d’equations differentielles du second ordre avec des termes iteratifs, on 
a utilise le theoreme de point fixe de Schauder. Et via celui de Banach et 
sous certaines conditions appropriees, on a etabli l’existence et l’unicite. 
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Chapitre 4. Existence et unicite des solutions periodiques d’une classe 
d’equations differentielles iteratives du second ordre 

4.1 Introduction 

Les equations differentielles a retard dependant de l’etat sont utilisees 
pour modeliser mathematiquement plusieurs phenomenes dans differentes 
disciplines scientifiques comme l’electrodynamique, l’automatique, la biolo- 
gie, l’epidemiologie et l’economie, etc. L’etude de ces equations a commence 
un peut avant avec le probleme preeminent de deux corps en electrodyna- 
mique et qui a reste dans les limbes jusqu’a l’apparition des travaux de Driver 
[22, 2E[ [22, 3D| dans les annees 60 du siecle dernier. 

L’equation differentielle fonctionnelle a retard iterative qui est une rela¬ 
tion entre une fonction, ses derivees et ses iterations, est un type particular 
des equations differentielles a retard dependant de l’etat. 

II semble que les travaux sur les equations avec termes iteratifs ont debute 
par un exemple simple mais tres important traite par Babage m en 1815 
ou il est interesse a trouver une fonction telle que sa n ieme iteration egale a 
la fonction elle meme. 

Dans le siecle suivant, ce type d’equations a attire l’attention de plusieurs 
chercheurs de l’epoque et une tres vaste litterature existe concernant l’etude 
de ce type d’equations par des techniques differentes. A titre d’exemples, on 
peut citer 

En 1968, P. Andrzej [5] a utilise la methode des approximations succes- 
sives de Picard , pour etudier l’equation 

x' (t) = f(t,x(t),x(x C t ))), 

ou 0 est un point de domaine extreme a gauche. 

En 1984, Par le principe de l’application contractante, Eder [32j a obtenu 
quelques resultats d’existence d’une solution monotone unique de l’equation 

x' (t) = x (x (t )), 
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avec x(t 0 ) = t 0 (^o e [-1, !])■ 

L’equation 

x' (t) = f(x{x (t))), 

avec x(a) = a, ou a est un point de domaine extreme d’un intervalle bien 
defini, a ete etudiee en 1990 par Wang [ST] et par M. Feckan j33] en 1993, 
ou x( 0 ) = 0 . 

Quatre ans plus tard, la meme equation a ete etudiee par W. Ge et Y. 
Mo [37] sous l’hypothese x(t 0 ) = x 0 dans un intervalle compact donne, ou 
les points extremes de l’intervalle sont deux points adjacents non nuls de /. 
Ce probleme a attire par la suite l’intension de Berinde [13] qui a applique 
la methode des operateurs non expansives dans son travail publie en 2010 . 
En 1995, Stanek [SS] a etudie les proprietes globales des solutions de l’equa- 
tion differentielle iterative suivante : 

x' (t) — x (t) + x {x (t)) 

et il a obtenu des resultats sur l’existence de solutions strictement monotones. 
En 1997, Si , Li et Cheng [711] ont considere l’equation : 

x' (t) = jW (t) ? 

et ils ont etabli des conditions suffisantes pour assurer l’existence de solutions 
analytiques. Apres une annee, Si et Wang [80] ont prouve l’existence des 
solutions regulieres de l’equation : 

x' (t) = Xrx (t) + A 2 x [2] (t) + ... + X n x M (t) + f (t) . 

Plusieurs chercheurs se sont penches ces dernieres annees sur l’etude de ce 
type d’equations en utilisant la technique du point fixe. Dans [ 02] , H. Y. 
Zhaoa et J. Liu ont applique les theoremes de point fixe de Banach et de 
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Krasnoselskii pour etablir l’existence, l’unicite et la stabilite de solutions 
periodiques de l’equation differentielle iterative du premier ordre suivante : 

= ci(t)x [1] (t) + c 2 (t)x^ ( t ) + ... + c n (t)x [n] (t) + F (t). 

Dans E5, H. Y. Zhaoa et M. Feckan ont utilise le theoreme de point fixe de 
Schauder pour montrer l’existence et la stabilite de solutions de l’equation 
suivante : 

k oo 

W = {x [m] (t)Y + G(t). 

m=11=1 

II faut signaler ici que la rrierne derniere equation a ete etudie par Liu et Si 
|5BJ en 2007 par la transformation de Schroder. 

Les travaux realises par la plupart des chercheurs jusqu’a maintenant res¬ 
tent concentres sur les equations du premier ordre, par contre on ne trouve 
que tres peu de travaux consacres a l’etude des equations du second ordre. 
Dans ce chapitre et via les theoremes de point fixe de Schauder et de Ba¬ 
nach, notre but majeur est d’etablir des conditions suffisantes pour prouver 
l’existence et l’unicite des solutions periodiques de l’equation differentielle 
iterative du second ordre suivante : 

+P(t) j t x(t) + q(t)x(t) = (t,x(t),x [2] (t),...,x [n] (t)) 

+h(t, x(t ), x [2] (t ),..., x [n] (t)) ,(4.1) 

ou p et q sont des fonctions positives continues, = t, = x(t), 

x^(t) = x (x(t )),..., = x 1 ™- 1 ^ (x(t)) et h, g : M n+1 —> M sont des 

fonctions continues par rapport a leurs arguments. 
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4.2 Preliminaires 


Pour tout T > 0 et L, M > 0, 

P t (L,M) = {x e Pt, \\x\\ < L, \x (t 2 ) — x (ti)| < M \t 2 — ti\ , Vti,t 2 e M} , 

est un sous ensemble ferme, borne et convexe de Pt- Les fonctions h(t, x\. x 2 -.... x n ) 
et g(t, Xi, x 2 , ..., x n ) sont supposees T—periodiques en t et globalement lip- 
schitziennes en X\, ...,x n i.e, 


h(t + T,x i, ...,x n ) = h(t,x i, ...,x n ) , g (t + T,x i, ...,x n ) = g (t,x i, ...,x n ) , 

(4.2) 

et qu’il existe n constantes positives k\. fc 2 ,..., k n et n constantes positives 
Ci, c 2 , c n telles que 

n 

\h(t,xi, ...,x n ) - h(t,y i, ...,|/ 2 )| < ||xi - yi\\ , (4.3) 

i =1 

et 

n 

I g (t,x !,..., x„) - g(t,y !, y 2 )| < _ 2/ill ■ ( 4 - 4 ) 

2=1 


Lemme 4.1 /3Hj/ Pour tout ip,7p E Pt (L, M ), 


72—1 


(/?["] — - 0^1 < M - 7 ||<yj — 011 , n = 1 , 2 ,.. 

j=0 


Lemme 4.2 Supposons que les hypotheses (3.2),(3.3), (4-2) et (3.7) soient 


verifiees. Si x E Pt (L, M ), alors x est une solution de (4-1) si et seulement 
si 


r t +T 


x (t) = I [E (t, s) — a(s)G(t, s)] g (s, x(s), x^(s), ..., x^(s)) ds (4.5) 


ft +T 


+ I G(t, s)h (s, x(s), x^(s),..., x^(s)) ds, 


ou 


E (; t , s) = 


exp(// 5 (u) du) 
exp(J Q r 5 (v) dv) — 1 


(4.6) 
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Preuve. 

Ligne directrice : En utilisant le Lemme 3.3 et certaines proprietes 
utiles de la fonction de Green citees dans le corollaire 3.1, on pent 
convertir l’equation (4.1) en une equation integrate. 


Supposons que x G Pt (L , M) est une solution de (4.1). Du Lemme 3.3 dans 
le chapitre precedent, il est facile de voir que 


rt +T 


x(t) = / G(t,s) 


d 


—g (s, x(s), x [2] (s),..., x [n] (s)) + h (t, x(t),x [2] (t), ...,x [n] 


rt +T 


G(t,s ) 


rt +T 


dt 


ds 


+ J G(t, s) [h (s, x(s), x^(s),..., x^(s))] ds. 
Une integration par parties nous donne 


rt +T 


G(t,s ) 


j t 9 (s,2;(s),x [2] (s),...,x w (s)) 


ds 


= [G(t,s)g (s,x(.s),a; l2l (.s),...,x w (s))] 


i t.+T 


rt +T / j 

—G(t, s) ) h (s, x(s), x^(s),..., x^(s)) ds. 


Comme 


II decoule du Corollaire |3.1 

rt +T 

/ G(t,s 


[G(t,s)g (s, x(s), x [2] (s),..., x [n] (s))] t = 0, 

que 

ds 


d 


—g (s,x(s),a; [2] (s),...,x [nl (s)) 


rt +T 


g (s, x(s), x^(s),..., x^(s)) [E (t, s ) — a(s)G(t , s)] ds. 


II s’ensuit que 


rt +T 


x (t) = I [E (t, s) — a(s)G(t, s)] g (s, x(s), x^(s),..., x^(s)) ds 


rt. +T 


+ 


G(t, s)h (t, x(t),x^ 2 \t ),..., x^(f)) ds. 


Ceci acheve la preuve du Lemme. 


(f)) ds 
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4.3 Existence de solutions periodiques 


Dans cette section, on va utiliser le theoreme de point fixe de Schau- 


der pour prouver l’existence de solutions periodiques de l’equation (4.1). Du 


Lemme 4.2 on definit un operateur S : Pt ( L , M ) —> Pf par 

rt +t 


{Sip) it) = 

+ 


[E (t, s) - a(s)G(t,s)]g (s, cp(s), cp l2] (s),..., (p M (s)) ds 


rt +T 


G(t,s)h (s,7>(s), v [21 (f>), ...,9 w W) ds, 


(4.7) 


et on va montrer qu’il est continu et compact et que Sip G Pr ( L,M ) pour 
tout ip G Pt (L, M) . Nous introduisons les notations suivantes : 


a = 


4 = 


T exp (Jq p(u) du) 


2 , £ = p v = exp 

1 


(e 1 - 1 r e — 1 v Jo 


exr 


p(/o a ( u ) du) — 1 exp( J q T 6 (it) du) — 1 


b (v) dv J , 

, Ai = max |a (t)| 


te[o,T] 


A 2 = max |6 (i)|, = max |/i (t, 0, 0,..., 0)|, p 2 = max |p (i, 0, 0,..., 0)| 

te[o,T] te[o,t] te[o,t] 

n j=i —1 n j=i—l 

( 1=Pl + Lj2k, (2 = P2 + Lj2 Ci E MJ - 

Z—1 j=0 z=l j=0 


(4.8) 


Lemme 4.3 Supposons que les conditions 1)3.2$ , 1)3.3 ), (3.1), (3.11) et (4-2)— (4-4) 


sont satisfaites. Alors Voperateur S : Pt ( L,M ) —> Pt donne par (4-1) est 
continu et compact. 


Preuve. 


Ligne directrice : La preuve est basee sur l’utilisation du 
theoreme d’Ascoli — Arzila afin de montrer la compacite 
de Pt (L, M) et partant la compacite de S decoule 
immediatement de sa continuite. 


85 
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Comme P T (L, M ) est un sous ensemble uniformement borne et equicontinu 
de l’espace des fonctions continues sur le compact [0, T ], alors le theoreme 
d’Ascoli— Arzila assure que Py ( L,M ) est un sous ensemble compact de cet 
espace. En outre, puisque tout operateur continu envoie tout ensemble com¬ 
pact en un ensemble compact, il nous suffit de montrer que S est continu. 
Pour tout ip, 9 e Pt (L, M ), on a 


\(S<P) (t) - (so) (01 

\E(t,s)\ g (s,<p(s),<p [2] (s),...,^ w (s)) -g (s,9(s),9 [2] (s),...,9 [n] (s] 


< 


ds 


ft +T 


+ 


+ 


a(s)| \G(t,s)\ g (s, 9(s), 9 [2] (s), ..., 0 [n] (s)) - g (s, 0(s), 9 l2] (s ),..., 0 w (s 


ft +T 




ds. 


Sous les conditions (4.3), (4.4) et d’apres le Corollaire 3.2, il vient 


\(Sp)(t)-(S9)(t)\<(^ + aX 1 )Tj2ci p [il -9 [i] PaT^h P [A - 9 [ 


i =1 


i=l 


En vertu du Lemme 4.1, il s’ensuit que 


j=i-l 


\(s«,)(t)-(se)(t)\<Tj2 ((/3 + aAi) Ci + aki) ^ M j \\p — 0||, 

3 =o 


i =1 


ce qui prouve que l’operateur S est continu. Done, S est compact. 


Lemme 4.4 Pour tout ti,t 2 £ 

rti+T 


\G(t 2 ,s)-G(t 1 ,s)\ds < \t 2 — ti| Te 2m g [T\ 2 u (2e 2m + l) + e m + l] 


ds 
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Preuve. Soit ti,t 2 G M, on a 
I G (t 2 , s) — G (ti, s)| 


< 


+ 


ft 2 ex p 

/“ 6 (v) dv + f* a (v) dv 

du — £ exp 

/“ 6 (u) cin + £ a (v) dv 

du 


exp(J 0 T a (w) <iw) — 1 

exp(J Q r b (u) du) — 1 



rt 2 +T 

Js eX P 

ft 2 b ( v )du + / U s+T a ( v ) dv 

du — £ 1+T exp 

/“ b (v) dv + /J +T a (v) dv 

du 


exp (£ a (w) du) — 1 

exp ( J q T b (u) du) — 1 



Mais 


exp 


' *2 


b ( v) dv+ a (v) dv 


-Jt2 


du — exp 
Jt\ 


b (v) dv+ a ( v) dv 


-Jti 


du 





S 2 b ^ dv x a(v)dv ) du- ( e /f “ fe(,;)ci,J x e £ dw 

t 2 V ' Jt 2 

2 'g/t” x e i: *(v)dv) du 

e fu a ( v ) dv ( e fZ b M*> - e ft u i b ^ d ^ du - f' 2 (e f *i b(v)dv x e fi a W dv ) du 



't2 


< 


I e s:«»)*> e f?Hv)dv ^ j : 2 b(v)dv _ e r tx &(«)*j du- J (e 1 ^ b{v)dv x e ^ s “W*) du 

rt 2 


^2 


= / t2 6 o)<^ _ p /ii &(«)<& 


e f* a (v)dv x e f^b(v)dv 


't2 


du + 


3 /*~ x e f:°Mdv du 


>t i 


< 

e ft 2 b ( v ) dv _ e ft\ b ( v ') dv 

f 

e If a(v)dv x e ff b(v)dv 

rt 2 

du + 

g If b(v)dv x e / 0 T a(v)dv 



Jo 


Jtx 



du 


< 


e ft 2 b ( v )dv _ e f t \ b(v)dv 


rt 2 


/ e /m dw + / 
'0 iti 


< Te 2m 

Comme 


e // 2 6 M*> _ e /ti 6 0)<*- 


+ e* m \t 2 -t 1 \. 


rti+T 

e ft 2 b ^ dv _ g/ij 

/*il “hX” 

ds= e ^)dv 

i _ g/A 2 

Iti 


Jti 

/ 

rT \ 



< T|| 6 ||e x p l-J 

b (v) dv ) \t 2 


ds 
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alors 


rti+T 

/‘ex pf 

It i 

Jti 


b (v) dv+ a (v) dv 


- ^2 


du — exp 
Jti 


b ( v ) dv+ a ( v ) dv 


-Jti 


du 


ds 



< re ||6|| exp _ / b (v) dv |t 2 — ti\ + Te \t 2 - h\ 


|t 2 -ti|Te 2m fr||6||exp 


D’une maniere similaire, on aura 



b (v) dv ) + 1 ) . 


(4.9) 


rt 2 +T 


e ft U 2 Kv)dv e f°+ T a (v)dv du _ 

f tl+T e ^ b{v)dv e^ +Ta ^ dv du 

S 

+ [ t2+T e f ** b(v)dv e^ +T a(v)dv du 
Jtl+T 

< I i+ e ItU 2 b(v)dv e^ +Ta{v)dv du 


-il +T 


git“ b ( v ) dv e Su +T a{v)dv d; 
tl+T . s+T 


U 


Jh b ^ dv e i: + <^ du 


rti+T 


Jt“ b ( v ') dv e fu +T a ( v ) dv d 


U 


+ 


rti+T 

\ 

'n+T 

rti+T 


fti+T 


^ b ^ dv e i: +T ^du 


e It 2 b ( v ) dv _ g/ti Kv)dv^ e /J+ T a(v)dv du 


+ 


r+i+T 


Jt 2 b ( v '> dv e fu +T a ^ dv du 


'il +T 


e ft 2 b(v)dv _ e f t \ b(v)dv\ e f°+ T a (v)dv g /“ b {v)dv du 


+ 


rti+T 


'tl+T 


,ft 2 b ( v ) dv e fu +Ta ( v ) dv d . 


U 


< 

+ 

< 


Jt 2 b ( v ) dv _ g/ti b ( v ) dv 


Jt 2 b ( v ) dv _ g/ij b ( v ) d ” ^ 1+ r2T 


[ 1+ e^ +Ta ^ v)dv e^ Kv)dv du 

+ 

/»£o —I— r [' 

f e /“ 2 b ^ dv e s: +T < v ) dv du 

Js 


Jt!+T 


e fo a H dv e fo b (v)dv 


du 


I e Io b ( v ) dv e fo T <*(v)dv 

ti+T 

Jt 2 b H dv _ Jt\ b (v)dv 


du 

2Te 4m + e 3m |t 2 -t 1 |. 
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Nous avons done 


rti+T 


't i 


r T e % b(v)dv e^ +T a ^ dv du - r e ; *i b{v)dv e C T <^du 

s J s 

T 

3m 


ds 


< 2T e H&ll exp — / b (v) dv \t 2 — t( + Te 6m \t 2 — h 


' o 


— \t 2 — t\\T [ 2 Te 


Am 


exp 


b (v) dv ) + e 


_3 m 


(4.10) 


II s’ensuit de ( |4.9[ ) et ( |4.10[ ) que 

rti+T 


|G ( t 2 , s ) — G (ti, s)| ds 


't i 


< |^2 — tl\T- 


2 Te 


4 m 


exp J 0 T b ( v ) dej + Te 3m + e 2m T ||b|| exp b (v) dv 


exr 


p(/o a ( u ) du) — 1 exp ( J Q r 6 (u) du) — 1 


= |*2 - *i| T 


T ||6|| exp / Q r b (v) dvj (2e 4m + e 2m ) + e 3m + e 2m 


exr 


p(J 0 a (u) du) — 1 exp(J Q r 6 (u) du) — 1 


< |t 2 - t x | Te 2r7 N [TA 2 z/ (2e 2m + l) + e m + l] . 


La preuve est achevee. 


Lemme 4.5 Supposons que les conditions (4-2)— ( 4 - 4 ) son t veriftees. Si 


T ((3 + aAi) Ci + Ta( 2 < L, 


(4.11) 


et 


[(2a + Te 2m q (T\ 2 v (2e 2m + l) + e m + l)) (AiC 2 + Ci) + (2/3 + TX 2 u) C 2 ] < M, 

(4.12) 

alors 

S (P T (L, M)) C P t {L,M). 


+ e 2m 
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Preuve. 

Ligne directrice : Pour assurer Vinclusion ci-dessus, il suffit done 
de montrer ces trois conditions : 

• S{(p)e P T , V<p e P T (L, M) 

• |(^) (t 2 )-(S0) (t 2 )\ <M \h-hl, Vt 1: t 2 e M, \/(p, 9 e P T (L, M) 

• \\S(<p)\\ <L, V<peP T (L,M). 


II est aise de voir que (Sip) (t + T) 
arrive a 


(Sip) ( t ). Grace au corollaire 3.2 


on 


r t +T 


|(<SV) (t )| < (P + aXi) 

pt +T 


+ Oi 


|c/(s,^(s),^ [2l (s),...,(^ [n] (s))|ds 

Jt 

|h (t,ip(t),ip m (t),...,ip [n] (t))\ds. 


Sous les conditions (4.2) —(4.4) et en vertu du Lemme 4.1, on obtient 


|h (s,^(s),^ [2] (s),...,^ [n] (s))| 

= |^(s,^(s),^ [2] (s),...,¥> [nI (s)) — h (s, 0, 0,..., 0) | 

+ |h ( s , 0 , 0 ,..., 0 )| 

n j—i— 1 n j=i—l 

< pi + Y ki y M1 imi <pi + LY k i Y MJ - £i> 

2=1 j= 0 2—1 j= 0 


et 


|g (s,¥»(s),^ [2] (s),...,¥> [n] (s))| 

= |c/(s,^(s),^ [2l (s),...,^ [n] (s)) — g (s, 0,0,..., 0) | 

+ 1 9 (s, 0, 0,..., 0)j 

n j=i— 1 n j=i—l 

<p 2 + Y Ci MJ imi <p2+ L Y Ci MJ - £ 2 - 

2=1 j =0 2=1 j =0 

Done 

l(^) (t)\<T((3 + a\ 1 )C 1 + TaC 2 . 
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De (4.11), il vient 


IISVII < L. 


Soit 1 1 , i 2 eR avec t\ < t 2 , on a 


I CSV) (t 2 ) - {Sip) (ti)| 

E{h,s)g (s,ip{s),p [2] {s),...,p [n] {s)) ds 


< 


n 2 

fii+T 


't i 


E(t u s)g (s,p{s),p [2] {.s),...,p [n] {s)) ds 


+ 


rt 2 +T 


't 2 
rti+T 


a{s)G{t 2 ,s)g (s,ip{s),p [2] {s), ...,p [n] {s)) ds 


+ 


Jti 

rt 2 +T 


't 2 


a{s)G{t u s)g (s, p(s), p l2, (s ),..., p M (s)) ds 
G(t 2 , s)h (s,p{s),p [2] {s),...,ip [n] (s)) ds 


fti+T 


'ti 


G(t i, s)h (s, <p(s), <p [2] (s), <p W ( s )) ds 


D’une part 


rti+T 


’t 2 


E{t 2: s)g (s,p(s),p [2] {s),...,p [n] {s)) ds 


r*ti+T 


'ti 


E {t u s) g (s, p{s), p [2] {s),p [n] {s)) ds 


< [ \E{t 2 ,s)\\g (s,<p(s),<p [2] (s),...,<p N (s))|cfc 

Jt 2 

rt 2 +T 

+ / \E {t 2l s)\\g (s,p(s),ip [2] (s),... 1 p [n] (s))\ds 

Jt x +T 

pti 

+ / \E{t 2 ,s) - Eit^s^lg (s,p(s),p [2] (s),...,p [n] {s))\ds 

Jtx 

< 2 \t 2 — t\\ (3(,2 + |^2 — t\\ T\ 2 vC t2 
= \h ~ ti \ (2/3 + T\ 2 u) ( 2 , 


(4.13) 
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et d’autre part 

ft2+T 


't2 


a(s)G(t 2 , s)g (s, if(s),^ [2] (s ),..., <^ w (s)) ds 


rh+T 


a{s)G{t 1 ,s)g (s, <p(s), </? [2] (s), <p [n] (s)) ds 


< 

+ 

+ 


f*t 1 


a(s)G(t 2 ,s)g (s, ip(s), ^ [2] (s),</? [n] (s)) ds 

Jt2 

/ a(s)G(t 2 ,s)g (s, tp(s), <^ [2] (s),<^ W (s)) ds 

>ti+T 
r^i+T 

/ a(s) [G(h,s) - G(t u s)}g (s,p(s),p [2] (s),...,ip [n] (s)) ds 

'ti 


Du Lemme 4.4 et les conditions (4.3) et (4.4), on aura 

rt2+T 

n 2 


E(t 2 ,s)g (s,p(s),ip [2] (s),...,p [n] (s)) ds 


r-ii+T 


E{ti 1 s)g (s,^(.s),^ [2] (.s),...,99 W (s)) ds 

J t\ 

< 2aXi \t 2 — ti| C 2 4" 1*2 — *i| ( 2 Te 2n \ (TX 2 v (2e 2m + l) + e m + l) 

= 1*2 - *i| AiC 2 (2a + Te 2m q (TX 2 u (2e 2m + l) + e m + l)) , (4.14) 


et 


rt 2 +T 


^2 


G(t 2 ,s)h (s,p(s),ip [2] (s),...,p [n] (s)) ds 


rti+T 


'ti 


G(t 1 ,s)h (s,ip(s),ip l2] (s),...,ip [n] (s)) ds 


< 

+ 

+ 




G(t 2 , s)h (s, p(s), <^ [2] (s),..., V5 N ( S )) ds 

^2 

rt2 +T 

/ G(t 2 , s)h (s,^(s),^ [2] (s),...,<^ w (s)) ds 

't!+T 
ft i~\~T 

/ (G(t 2 ,s) - G(t 1 ,s))h(s,ip(s),^ l2] (s),...,^ [n] (s)) ds 

'ti 
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Grace au Lemme 4.4 et les conditions (4.5) —(4.6), on arrive a 

r t2+T 


>t 2 


G(t 2 ,s)h (s,ip(s),ip [2] (s),...,ip [n] (s)) ds 


rti+T 


'ti 


G(t u s)h (s,^(s),<p [2] (s),...,<p w (s)) ds 


< 1 t 2 - til Cl (2a + Te 2m q (T\ 2 v (2e 2m + l) + e m + l)) 


(4.15) 


Ainsi, il s’ensuit de (4.13), (4.14) et (4.15) que 


I (S<p) (t 2 ) - (Sip) (ti)| 

^ |^2 — £i| (( 2 a + Te 2m q (TX 2 U (2e 2m + l) + & m + l)) (A1C2 A Ci) A ( 2/5 + T^v) C2) 


Moyennant la condition (4.12), on obtient 


|(<svo (t 2 ) - (Sip) (ti)| < m 1 1 2 - ti 

Ce qui acheve la preuve. ■ 

Void maintenant notre resultat principal : 


Theoreme 4.1 Supposons que (3.2), (3.3), (3.2), (3.12), (4-2)— (4-4' (4-H) 


et (4-12) sont veriftees, alors l’equation (4-1) admet au moins une solution 


x € Pt (L, M). 

Preuve. 

Ligne directrice : Le resultat decoule immediatement 
des Lemmes 14.31 et 14.51 


Du lemme 4.2[ on voit, que l’equation (|4.1|) admet au moins une solution 

x 


dans Pt ( L , M) si et seulement si l’operateur S defini par ( 4.7[ ) admet au 
moins un point fixe. Pour conclure, il suffit d’appliquer les Lemmes |4.3| et |4.5| 
ce qui prouve que les hypotheses du theoreme de Schauder sont satisfaites. 
Par consequent, il existe au moins un point fixe x G Pt (L, M ) de l’operateur 
S. D’ou l’existence d’au moins une solution de l’equation (4. 1[) . ■ 
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4.4 Existence et unicite des solutions perio¬ 
diques 

Nous donnons pour clore ce chapitre un theoreme qui assure l’existence 


et l’unicite des solutions de F equation (4.1). 


Theoreme 4.2 Supposons (3.2), (3.3), (3.17), (4-2)— (4-4) sont satisfaites. 
Si 

n j=i— 1 

” " (4-16) 


T ^ ((/5 + aAi) Cj + aki) ^ M j < 1, 

i=o 


i —1 


alors l’equation (4-1) admet une solution unique x E Pt ( L,M ). 


Preuve. Soit <p, 6 G Pt (L, M ). Par une demonstration analogue a celle du 


Lemme 4.3, on obtient 


j=i -1 


((/3 + aAi) Ci + aki ) AF J ||<p — 0||, 

l=o 


2=1 


ce qui entraine que moyennant Fhypothese (4.16), S est une contraction. 


Sous la condition citee ci dessus, et en vertu du theoreme de point fixe de Ba¬ 
nach, S admet un point fixe unique dans Pt ( L , M), ce qui prouve l’existence 


d’une solution unique de Fequation (4.1) 
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4.5 Recapitulatif 

Nous avons aborde dans ce chapitre quelques questions liees a l’existence, 
l’unicite et la periodicite des solutions d’une classe d’equations differentielles 
fonctionnelles iteratives du second ordre. Les resultats sont obtenus en se 
basant essentiellement sur une methode combinant le theoreme de point fixe 
de Schauder et les proprietes de la meme fonction de Green du chapitre 
precedent. 

Nous montrons aussi que, sous certaines conditions et en vertu du theo¬ 
reme de point fixe de Banach , la meme equation differentielle iterative pos- 
sede une solution unique. 
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Conclusion generale et Perspectives 


L’interet de la theorie de point fixe pour l’etude des equations fonctionnelles 
non lineaires a attire mon attention et c’est pour cette raison que je l’ai choisi 
comme sujet de ma these de Doctorat ou nous sommes interesses a etablir 
des conditions sufhsantes pour assurer P existence, l’unicite, la periodicite et 
la positivite des solutions de deux classes d’equations different idles a retard 
variable du second ordre ou la premiere classe, est une classe d’equations de 
type neutre tandis que la deuxieme est d’equations de type iteratif. 

Au cours de ce travail, nous avons presente un etat de Part sur les equa¬ 
tions fonctionnelles a retard et nous avons egalement mis en evidence la 
grande efhcacite de la technique qui combine quelques proprietes utiles des 
fonctions de Green de certaines EDO lineaires avec les theoremes de point 
fixe. 

Cette technique puissante repose sur Pinversion de l’equation consideree 
en une equation integrate avant d’appliquer le theoreme de point fixe de Kras- 
noselskii ou celui de Schauder combine avec quelques proprietes d’une fonc- 
tion de Green d’une equation differentielle lineaire du second ordre. En outre, 
nous avons applique le principe de l’application contract-ante de Banach sur 
l’equation integrate inversee et sous certaines contraintes determinees, on a 
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Conclusion et perspectives 


demontre l’existence et l’unicite de la solution periodique de l’equation dif¬ 
ferent ielle de type iter at if. 

Ce type d’equations est a la mode aujourd’hui et un changement soit 
dans l’espace de depart ou soit le type du retard ouvre la voie vers pleusieurs 
autres pistes interessantes ayant trait au points suivants : 

— II apparait interessant de demontrer la stability de ces solutions. 

— II serait aussi tres utile de pouvoir etendre l’etude aux equations diffe¬ 
rent ielles a retard d’ordre superieur a deux, equations differentielles a retard 
distribue ou equations fractionnaires a retard. 

— II semble importent d’etudier l’existence de solutions presque perio- 
diques ou anti periodiques. 
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Annexe 


I. Resolution des EDR par MATLAB 

Les principales fonctions de MATLAB qui permettent la resolution des equa¬ 
tions differentielles a retards sont dde23, ddesd et ddensd. 

1. Resolution d’une equation differentielle a retard constant par 
dde23 

Le solveur dde23 de Matlab, base sur les methodes de Runge-Kutta 
explicite est developpe par S. Thompson et L. Shampine pour resoudre 
les equations different ielles avec des retards constants. 

Syntaxe de base 

sol = dde23 ( ddefile, lags, history, tspan) 

sol = dde23 ( ddefile, lags, history, tspan, options), 

ou les arguments d’entree sont definis par : 
ddefile : fonction au second membre de l’EDR. 
lags : vecteur des retards constants positives. 

history : fonction historique, c’est un vecteur colonne donnant les 
valeurs initiales des solutions. 

tspan : intervalle d’integration (temps ou on veut calculer la solution). 
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2. Resolution d’une equation differentielle a retards variables de¬ 
pendant continument du temps ou de l’etat par ddesd 

La fonction ddesd permet la resolution des equations differentielles avec 
des retards variables. 

Syntaxe de base 

sol = ddesd (ddefun, delays, history, tspan) 

sol = ddesd (ddefun, delays, history, tspan, options), 

ou les arguments d’entree sont definis par : 
ddefun : fonction creee qui evalue le second membre de l’EDR 
delays : vecteur colonne qui represente les retards qui peuvent- de- 
pendre du temps ou de l’etat. 

history : fonction historique, c’est un vecteur colonne donnant les va- 
leurs initiales des solutions. 

tspan : intervalle d’integration (temps ou on veut calculer la solution). 

3. Resolution d’une equation differentielle a retard de type neutre 
par ddensd 

La fonction ddensd permet la resolution des equations different ielles de 
type neutre avec des retards independants de l’etat. 

Syntaxe de base 

sol = ddensd (ddefun, dely, delyp, history, tspan) 

sol = ddensd (ddefun, dely, delyp, history, tspan, options). 

ou les arguments d’entree sont definis par : 

ddefun : fonction creee qui evalue le second membre de l’EDR 

dely : vecteur colonne qui represente les retards. 

delyp : vecteur colonne qui represente les retards dans les derivees. 
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history : fonction historique, c’est un vecteur colonne donnant les 
valeurs initiales des solutions. 

tspan : intervalle d’integration (temps ou on veut calculer la solution). 

Un tutoriel sur la resolution des equations differentielles a retard via 
MATLAB par Shampine et Thompson peut etre trouve sur le site 

Webhttp ://www.runet.edu/thompson/webddes/ 

II. Resolution des EDR par Maple 

La premiere procedure est appelee fhmpas, repose sur une resolution 
exacte des EDO successives. Elle devient lente des que le nombre d’iterations 
devient grand et doit etre alors remplacee par des procedures de resolution 
approchee, appelees fhmpasn et fhmpasndp , utilisant les ressources de l’option 
numeric dans dsolve. 
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